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D , Abstract 
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We develop a systematic and efficient method of counting single-trace and 
multi-trace BPS operators with two supercharges, for world-volume gauge theories 
of N D-brane probes for both N — > oo and finite N. The techniques are applicable 
to generic singularities, orbifold, toric, non-toric, complete intersections, et cetera, 
even to geometries whose precise field theory duals are not yet known. The so- 
called "Plethystic Exponential" provides a simple bridge between (1) the defining 
equation of the Calabi-Yau, (2) the generating function of single-trace BPS oper- 
ators and (3) the generating function of multi-trace operators. Mathematically, 
fascinating and intricate inter-relations between gauge theory, algebraic geometry, 
combinatorics and number theory exhibit themselves in the form of plethystics 
and syzygies. 
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1 Introduction 

The study of BPS states in a quantum field theory is of unquestionable importance. 
The purpose of this note is to discuss the set of all mesonic BPS gauge invariant oper- 
ators (GIO) with two supercharges which appear in the chiral ring of a generic Af = 1 
super symmetric gauge theory that lives on a D3-brane which probes a singular Calabi 
Yau (CY) manifold 1 . For arbitrary singularities, finding the gauge theory living on the 
D3-brane is intricate. The simplest class is the orbifolds, the study of which began with 
[2, 3, 4, 5, 6, 7, 8]. The next simplest class is the toric singularities, the investigation of 
which was initiated by [8, 9, 10, 11]. Interesting duality structures of these theories have 
been expounded in [12, 13, 14, 15]. It is recently realised that the toric theories are, 
in fact, best described using a bi-partite periodic tiling of the two dimensional plane, a 



so-called "dimer model" [16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28] (also cf. [29, 30] 



for recent mathematical treatments). 

When the manifold is non-orbifold and non-toric there is no current systematic way 
of describing the gauge theory even though some examples exist in the literature. For 
example, the higher del Pezzo series [31] and certain deformations of toric singularities 
[27] have been addressed. In this paper we will see how one can describe the single- and 
multi-trace operators in terms of generating functions which can be computed for both 
1 Somc preliminary results were annouced in [1]. 
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toric and non-toric manifolds. In fact, the computations we will see can calculate the 
generating functions even for cases in which the gauge theory is not precisely known - 
either the superpotential is missing or even the quiver itself is not known. 

The discussion on GIO's in the chiral ring can be divided into few parts as follows. 
Given a gauge theory description of the theory on the D-brane, there are several problems 
of interest: 

Global U(l) Charges: One would like first to identify the set of global U(l) charges 
of this theory. One charge out of this set is singled out to be the R-charge and the other 
charges can be generically called global non-R charges. The most useful way of thinking 
about these charges is by introducing the holographically dual gravity description. A set 
of D3-branes on a singular conical CY is holographically dual to an AdS 5 x Y 5 background 
where Y 5 is a Sasaki Einstein (SE) manifold (cf. [32] and references therein). The global 
charges of the gauge theory are dual to gauge fields in AdS 5 . These gauge fields can be 
divided into two sets - one set originates from the Type IIB metric, those are typically 
referred to as the isometry of the SE manifold, and the other set comes from the Type 
IIB 4-form. The R-charge is always part of the isometry group of the SE manifold. 
The traditional name which was given to the charges coming from the metric are flavor 
charges and those which come from the 4-form are called baryonic charges. 

The isometry group of the SE manifold has a maximum rank of 3, in which case the 
SE manifold and its CY cone are called toric; the minimum rank is 1, in which case the 
corresponding U(l) charge is precisely the R-charge. The number of baryonic charges 
is in principle unbounded and is given by the third homology of the SE manifold. Most 
cases which were studied in the literature have one baryonic symmetry, the prototypical 
example being the conifold [33]. Currently there are extensive studies of cases with 
more than one baryonic charge, the simplest being the Suspended Pinch Point (SPP) 
[9, 18, 34], as well as the more complicated X p ' g family [50]. 

Counting Gauge Invariant Operators: Given the set of U(l) symmetries, R, 
flavor, baryonic, etc., say, n of them, each gauge invariant operator in the chiral ring 
carries a set of charges under these symmetries. We will assign a generic complex variable 
ti,i = 1 . . .n to each such charge and define a function f{{ti}) to be the generating 
function of all these operators. This function / has, by definition, an expansion in 
terms of monomials in {tj} such that the coefficient, c fel) ... )fen of t^ 1 • • -t^ n is integer and 
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counts the number of operators of charges (k±, . . . , k n ), 

MU})= E ^,.,kJ k l l ---t k n ■ (1-1) 

h,—,ik 

Our goal is to compute such functions for a multitude of cases. 

Our ultimate wish is that for any CY manifold we would like to know 

1. The set of single-trace BPS operators, the generating function is denoted by /; 

2. The set of multi-trace BPS operators, the generating function is denoted by g; 

3. For N D3-branes at the singular CY we would like to know the dependence on N. 
Namely, we would like to know how many independent single-trace and multi-trace 
operators are there in the chiral ring for a given set of charges. For a finite N this 
turns out to be a much more difficult task since there are matrix relations for a 
finite size matrix that need to be taken into account. Nevertheless, we propose a 
nice solution to this as well; the generating functions in this case will be denoted 
as Jn for single-trace and qn for multi-trace. 

As will be discussed in detail in later sections there is an important function which 
beautifully relates the single-trace generating function and the multi-trace. Namely, 

can be simply computed from f N using the so-called "Plethystic Exponential." This 
function has been used in physics several times in the past and we believe it should go 
into the literature more often as it plays a crucial role in counting problems such as 
the one dealt with here 2 . Conversely one can use the so-called "Plethystic Logarithm" 
which is the inverse function to the plethystic exponential and computes fx from g^. 
The ability to switch between / and g will turn out to be a key tool in analyzing the 
theories we are interested in and to reveal new pieces of information which were either 
previously unknown or alternatively not well discussed. 

Having presented a host of functions and concepts, it would be most expedient to 
exemplify them in a context with which the readers are well-acquainted. We shall do 
so for the famous D3-brane theory on C 3 in the next section. Having whetted the 
readers' appetites, the plan for the remainder of the paper is as follows. We begin 
with the large N limit and present the solution to questions (1) and (2) above. In §3, 
we show how to construct /, the generating function for single-trace GIO's. This is 
2 A. H. would like to thank Marcos Marino for demonstrating the properties of this function [53]. 
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a Hilbert-Poincare counting problem. We exemplify with orbifolds, toric varieties and 
the del Pezzo family. We take an interlude in §4 and examine this counting problem 
using the graphical perspectives of dimers. Then, in §5, we construct the generating 
function g, which count the multi-trace GIO's. The relationship between / and g will 
turn out to be a plethystic one. In due course, we will show how plethystics actually 
encode not only the GIO's counting, but also the defining equation of the singularity. 
Interesting partition identities as well as syzygies in graded polynomial rings emerge. 
Having constructed the generating functions, we then calculate the asymptotic behaviour 
thereof in §6, using results from combinatorics and analytic number theory. Finally, we 
use the above formalism to address the more difficult problem of finite iV in §7 and show 
how plethstics again solves the counting problem and how they encode the geometry. 
We conclude with perspectives in §8. 

2 C 3 : An Illustrat ive Example 

As promised in the introduction, we begin with a familiar example to illustrate the 
various generating functions. Here, the computation can be done without using the 
more general techniques which will follow in the rest of the paper. This example is of 
course for the archetypal example of the AdS/CFT correspondence, the case in which 
the CY manifold is trivially C 3 and its associated SE manifold, S* 5 [35]. There are no 
baryonic charges in this case since the third homology of S* 5 is trivial and the isometry 
group is SU(4) with rank 3, meaning that this CY manifold is actually toric and the 
number of U(l) charges is 3. We can thus define 3 corresponding variables, £i, £2, ^3, 
which will then measure these three U(l) charges in their powers, as explained above. 
The gauge theory is the M = 4 gauge theory with U (N) gauge group which in Af = 1 
language has 3 adjoint chiral multiplets which we will denote as x, y and z. Being toric, 
this CY manifold admits a description in terms of periodic bi-partite tilings of the two 
dimensional plane and in fact is given by the simplest of them all - tilings by regular 
hexagons [17]. 

We are interested in operators in the chiral ring and therefore we need to impose 
the F-term relations coming from the superpotential W = Tr(x[y,z]). The F-terms 
hence take a particularly simple form: [x,y] = [y,z] = [z,x] = 0, i.e., all chiral adjoint 
fields commute. The generic single-trace GIO in the chiral ring will then take the form 
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of Tr(x l y^ z k ). It is then natural to assign t\ as counting the number of x fields, 
the number of y fields and t 3 , the number of z fields in a GIO. There will therefore 
be a corresponding monomial t\t J 2 t k for each gauge invariant of charges and k, 
respectively. In fact, there will be precisely one for each triple of charges, provided each 
of and k are non-negative. Putting all of this together, we find that the generating 
function / takes the form 



00 00 00 



fitu^t^c 3 ) = YYYA4A = 7 w 1 w 7 • (2-i) 

' 123 (i-t 1 )(i-t 2 K J 

i=0 j=0 k=0 v in ZJK 61 

To be more precise, in the above form we did not take into account any relations that 
a finite matrix should satisfy, therefore, as mentioned earlier, this result is strictly valid 
for the case of N = 00. Therefore, using the notation introduced above we should write 

/oo(«i,« 2 ,*3;C 3 ) = (1 _ ti)(1 ^ t2)(1 _ t3) - (2.2) 

A General Feature for Toric CY: Note that in eq. (2.2) the coefficients c^fe 
appearing in the general expansion 

foo(tl, t2, fa) = Cjjkt\t^t\ (2.3) 
ijk 

are all equal to either 1 or 0. This means that for a given set of charges, k, there 
is either one operator carrying these charges or not, but there can not be more than 
one. Indeed this is a generic feature which is obeyed for every toric singular CY. More 
explicitly there is a one-to-one correspondence between single-trace GIOs and integer 
lattice points in the dual cone of toric diagram [10, 27, 37]. This property is reminiscent 
of some kind of a fermionic degree of freedom that carries this set of charges. In contrast, 
for the non-toric case, it is shown in [27] that there are, in general, multiple-to-one 
mappings between single-trace GIOs and given charges. The reason is clear. In the 
toric case, we have two extra U(l) flavor symmetries besides the R-symmetry, which 
is big enough to distinguish finely, while for non-toric case we do not have these extra 
symmetries. 

Let us look at the set of charges i,j,k for which are not zero. They form a 
sub-lattice of the three dimensional lattice which has the form of a cone. Indeed, this 
sub-lattice is the so-called "positive octant" for which i > 0, j > 0, k > 0. This feature 
of a cone structure will also be general for every CY manifold, the form of this cone is 
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interesting and will be discussed in detail in §4. One can think of the function as 
a theta-function over the lattice points of the cone and is a characteristic function of 
this cone. It is worth to notice that from results in [27], it seems that there is a lattice 
structure for both toric and non-toric cases. The difference is that for the toric case the 
lattice is 3-dimensional while for non-toric case the dimension is lower. 

If on the other hand we are interested in counting the number of BPS operators 
which carry a given fixed scaling dimension, say Tr(x l y^z k ) of dimension i+j + k = 
we need to set t-i — i 2 — t 3 — t in eq. (2.2) and get the generating function for all 
operators. In other words, we have to forget the other two U(l) flavor symmetries and 
use the fact that all variables x,y,z have same R-charge |. Hence, 

v ' m=0 v 7 

and the number of GIO's of given R-charge R = |m is ( m ^ 2 ), corresponding to the 
completely symmetric rank m representation of SU (3) that acts on x, y, and z in the 
fundamental representation. 

Single- Trace and Multi- Trace at N — > oo: Having studied /oo, let us now look 
at the function fi, generating the single-trace operators for the case of one D3-brane on 
C 3 . Clearly, the adjoint fields x,y, and z are complex variables and not matrices and 
therefore any product of two or more of these matrices is a multi-trace operator. As a 
result, there are only 4 single-trace operators in this case: the identity operator, x,y, 
and z. We can therefore use their representation in terms of t^i = 1,2,3, sum them 
and write: 

fi(ti,t 2: t 3 ) = l + t 1 + t 2 + t 3 . (2.5) 

Next, we notice an interesting relation between g 1 and /oo- Let us look at the set 
of operators of the form Tr(x l y j z k ) for the case in which the number of D3-branes is 
N — > oo. Each such operator is represented by the monomial t\t J 2t 3 and can be thought 
of as a multi-trace operator for the case of the number of D3-branes being N — 1. This 
implies that g 1 , the generating function for multi-trace operators for one D3-brane is 
equal to /oo, the generating function for single-trace operators for infinitely many D3 
branes, 

9i = /oo ■ (2.6) 
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Can we now find some functional dependence between f\ and g{l Combining ex- 
pressions eq. (2.6), eq. (2.5) and eq. (2.2), we have 

9i(ti,t 2 ,t 3 ) = _ j _ = exp[- log(l - h) - log(l - t 2 ) - log(l - h)} 
= expff^illLe^ff; «WH). (2 . 7 ) 



r=l ' x r=l 



The last relation 



g 1 {t l MM)=^(f2 mA f ) - ) =/oo(*l,*2,*3) (2-8) 

turns out to be a key relation and is satisfied for any CY manifold, toric or otherwise. 
The function g\ is then said to be the Plethystic Exponential of j\. This relation in 
fact generalizes to any N and we find that is the plethystic exponential of f^. We 
will discuss this extensively in §5 and §7. 

We are now ready to write down the expression for the generating function goo of 
multi-trace BPS GIO's in the chiral ring in the M = 4 theory, corresponding to iV — > oo 
D3-branes on C 3 . It is again, the plethystic exponential, this time of in eq. (2.2): 

*.(«.. t,, h) = ex P ( f; * - * ) = exp ( f; ~ i) . (2 . 9) 

Note that has an expansion 

9oo(ti,h,t3) — djj k t\t{t^ (2-10) 

where the coefficients are non-zero precisely when the coefficients of are 
non-zero. However, while can be at most 1, d^ has a very fast growth and in fact 
grows exponentially. It is therefore a problem of interest to find what is the large charge 
behavior of d^. We see that the multiplicity of BPS states for fixed R charge, R = 
is 

(OO 1 ^ \ oo 

— i=J2 d ^ k - ( 2 - n ) 
r=l ' k=0 

We will present in §6 detailed discussions of how to obtain d k for large k. 

Single- Trace and Multi- Trace at Finite N: For finite N, the situation is in gen- 
eral much more involved. Nevertheless, the multi-trace result g^ can be obtained from 
In by plethystics. We will present the systematic treatment for arbitrary singularities 
in §7. 
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3 Counting Gauge Invariants: Poincare Series and 
Single- Trace 

Having stated our problem and enticed the reader with the example of C 3 , we are now 
ready to attack the general CY singularity. Our strategy will be to first examine the 
simpler case of N — > oo and then the more involved case of finite N. 

Beginning with the large N situation, we first find the generating function / for 
the single-trace GIO's. Then, in §5, we will show how the plethystic exponential (PE), 
extracts g, the generating function for the multi-trace GIO's, from /. Indeed, because 
the multi-trace GIO's are composed of products of the single-trace ones, PE is expected 
to be a version of counting integer-partitions. We would like to emphasize that the 
counting automatically encodes more than merely the matter content, but, furtively, 
the superpotential as well. In other words, we will be concerned with a true counting 
of the GIO's with the F-term constraints automatically built in. We will check in all 
examples below that this is indeed so by showing that the moduli space is explicitly the 
CY 3-fold, as is required in D-brane probe theories. 

How, then, do we compute / given the geometrical data of the CY? It turns out 
that we could appeal to some known methods in mathematics. In projective algebraic 
geometry, an important problem is to count the number of generators of graded pieces 
of polynomial rings, the generating functions of this type are called Hilbert-Poincare 
series. 

We shall borrow this terminology and refer to the function / for the single-trace 
GIO's as the Poincare series for the associated M = 1 gauge theory; it shall soon be 
seen that this appropriation is a conducive one. In this section, we proceed stepwise 
along the various known classes of CY singularities which the D3-brane can probe. We 
start with orbifolds and see the Poincare series in the mathematical sense is precisely 
what is needed. Next, we address toric CY singularities; here, using the techniques of 
{Pi <?)-webs and 2-dimensional tilings (dimers), we construct / from the toric diagram. 
Then, we proceed to the del Pezzo family of singularities. 

3.1 Orbifolds and Molien Series 

Given a finite group, it is a classical problem to find the generators of the ring of 
polynomial invariants under the group action. The theory matured under E. Nother 
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and T. Molien (cf. e.g. [38]). In our quiver gauge theory, the single-trace GIO's are 
polynomial combinations of fields which are invariant under the group action. Because 
we are assuming large N, no extra relations arise beside these from the F-terms, and 
the problem of computing / reduces to simply counting the number of algebraically 
independent polynomials one could construct of degree n that are invariant under the 
group. The problem is a mathematical one and was solved by Molien; the Poincare 
series is named Molien series in his honour. 

Let us be concrete and specialise to the orbifolds of our concern, viz., 3-dimensional 
CY orbifolds C 3 /G, with G a discrete finite subgroup of SU(3). Such singularities 
were first classified by [39] and the D-brane quiver theories, constructed in [6]. Let G 
act on the coordinates (x,y,z) of C 3 . Then, the question is: how many algebraically 3 
independent polynomials are there of total degree n in (x,y,z). The Molien series is 
given by 



where the determinant is taken over the 3x3 matrix representation of the group elements. 
Upon series expansion, the coefficients bi give the number of independent polynomials 
in degree %. Hence, the / we seek is simply M(t; G). 

We can remark one thing immediately. In eq. (3.1) there is only one variable t instead 
of (x,y,z) in our example C 3 . The reason is that for orbifold theories which descend 
from the M = 4 parent every elementary field has R-charge 2/3. The replacement 
x,y,z — > t tells us that eq. (3.1) counts the single-trace GIO for given R-charge. Indeed, 
as a first check, take G = I, the trivial group. We immediately find that 



which agrees with eq. (2.2) for the C 3 theory if one set U = t. Thus, the Molien series 

counts invariants of total degree in x,y,z whereas eq. (2.2) counts the degree of the 

three variables individually. In the next subsection, we shall refine the Molien series by 

straight-forwardly generalising the dummy variable t to a triple (£1,2,3)- 

Emboldened by this check, let us go on to a non-trivial example, the binary dihedral 

group D 4 of 8 elements. This is a subgroup of SU(2) C SU(3) and is a member of 

the ADE-series of CY two-fold (K3) singularities (cf. [40]). We can think of this as a 

3 In fact linearly independent, because any polynomial relation would change the total degree. Find- 
ing the polynomial relations is a important one and will be subsequently addressed. 




(3.1) 



M(t; I) 



1 1 



1 + 3 1 + 6 1 2 + 10 t 3 + 15 t A + 21 1 5 + 0(t 6 ) , (3.2) 



det(I-tl) (1-t) 3 
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C 3 orbifold with one coordinate, say z, held fixed. The gauge theory is the well-known 
Af = 2 D-type quiver (q.v. [3, 4, 5, 6]). 

This group is generated (we use the standard notation that (xi, . . . , x n ) is the finite 
group generated by the list of matrices 



This dictates that there are two invariants at degree 4, one at degree 6, etc. 

Now, one can actually determine the invariants explicitly, which gives us another 
check. First, an important theorem due to Nother states that (cf. e.g. [38]): 

THEOREM 3.1. The polynomial ring of invariants is finitely generated and the degree 
of the generators is bounded by \G\. 

Therefore, though the Molien series is infinite, with increasingly more invariants 
arising at successive degree with them being linearly independent at each total degree, 
there will be highly non-trivial algebraic relations amongst the ones at different degree. 
The power of Theorem 3.1 is that one needs to find invariants at most up to degree 
equal to the order of the group; all higher degree invariants are polynomials in these 
basic ones. 

Hence, we only need to find a finite number of invariants, which can be determined 
explicitly due to an averaging technique of O. Reynolds (cf. e.g. [38]). Given any poly- 
nomial F(x), one can define the so-called Reynolds operator 



Then, the polynomial Rq[F(x)] is invariant by construction. We can then list all mono- 
mials of a given degree, apply eq. (3.5) to each and obtain the invariants at the said 
degree; the number thereof should agree with what eq. (3.1) predicts. 




(3.3) 



acting on (x,y) G C 2 . We can readily compute the Molien series to be 





(3.5) 
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Applying the above discussion to our example of D4, we obtain the following invariant 
polynomials for the first few degrees: 



degree 


invariant polynomials 


4 


*V, \{x A + y A ) 


6 


\xy(x A - y A ) 


8 


xY: ±*V(* 4 + 2/ 4 ), U x8 + V 8 ) 



(3.6) 



We remark that there are no invariants of lower degree (except trivially the identity) 
and that the number of independent invariants indeed agree with the series expansion 
of eq. (3.4). Immediately, one sees some trivial relations such as x A y A = (x 2 y 2 ) 2 . Using 
Grobner basis algorithms [42], one can show that the above ring of 6 invariants can be 
further reduced to 3. In other words, the ring of invariant polynomials, C[x, y] D4 , is 
generated by 3 so-called primitive ones: 

v = ^ 4 + y A ), w = x 2 y 2 , u = ^xy(x A - y 4 ) . (3.7) 

Finding relations among these polynomials is known as the syzygy problem and is, 
again, a classical problem dating to at least Hilbert. The modern solution is, as above, 
to use Grobner bases. The reader is referred to [43] for a pedagogical application of 
syzygies and Grobner basis to M = 1 gauge theories and to [44] in the context of moduli 
stabilisation. We will return to syzygies later in the paper. For the present example, we 
find the relation 

v 2 w - w 3 = u 2 . (3.8) 
This is a comforting result. Indeed, invariant theory tells us that 

The defining equation of an orbifold is the syzygy of the primitive invariants. 

We recognise eq. (3.8) as precisely the defining equation [40] for the afline variety C 2 /L>4. 
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3.1.1 ADE-Series 



For completeness, let us compute (making extensive use of [42, 45]) the Molien series 
for the discrete subgroups of 577(2). We find that 



G c £77(2) 



G 



Generators 



Equation 



Molien M(t; G) 



A 



n-l 



n 



UJ r . 







u- 1 



uv = w 



:i+f 



D n 



+2 



4n 



u 2r 




\ i 



y 2n 



% 



u 2 + f 2 w = w n+1 



(l+gn+2j 
'l-t 4 )(l-t^ 



l-t 4 +t 8 
l_ t 4_ t 6 +t lU 



24 



(5,T> 



■u 2 + f 3 + w 4 = 



£ 7 



48 



(5, 10 



■u 2 + + ww 3 = 



l-t 6 +t 12 
l_ t 6_ t 8 +t 14 



Ea 



120 



u 2 + v 3 + w 5 = 



l +t 2_ t 6_ t S_ t W +t U +t 16 



l +t 2_ t 6_ t 8_ t 10_ t 12 +t 16 +t 



18 



(3.9) 



where we have defined u„ := e » and 



.S':=i[ + * T 
1+i -1-i/ 

i f 1 + i , 



V2 



1-i. 



i 
.0 



4 4 4 



l-v 7 ! _ • l+v / 5 N 
4 4 4 

2 

2 



(3.10) 



We have also used the identity 



71—1 j 



oo / oo 



E E ^' +m ™w =«EE * 2m+n " + E < 



2m—n(3 



j=0 m=0 

// /' 1 



m=0 \/3=0 



1-t 2 V l - t n + t 



(3.11) 



3.1.2 Valentiner: A Non-Abelian SU(3) Example 

Having warmed up with the 2-dimensional CY orbifolds, we are ready to study the 
proper subgroups of SU(3) [6]. The simplest, most well-known, non-trivial, non-Abelian 
discrete subgroup of SU (3) is perhaps the Valentiner group, otherwise known as A(3-3 2 ) 
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(or, sometimes known as the Heisenberg group for 3 elements, as recently studied in [41]), 
defined as 
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(3.12) 



Let us investigate this group in some detail; we shall return to this group later in the 
paper. The Molien series is readily computed to be 

-1 + t 3 - t 6 



M(f; A(27)) 



l + 2r + 4t b + 7t 9 + llt i2 + 16t 15 + 22t i8 + ... . (3.13) 



(-1+t 3 ) 3 

To find the defining equation (syzygies), Theorem 3.1 tells us that we need only 
go up to degree 27 here, a total of 174 invariants, of degrees 0, 3, 6, ... , 24, 27. Using 
Grobner techniques [42], we find that there are only 4 nontrivial generators for these 
174 polynomials (we have scaled the expressions by an over-all 3): 

[m = 3xyz, n = x 3 + y 3 + z 3 , p = x 6 + y 6 + z 6 , q = x 3 y 6 + x 6 z 3 + y 3 z 6 } . (3.14) 

We then find a single relation in C[m, n,p, q\: 

8m 6 + m 3 (-48 n 3 + 72np + 72 g) + 81 (^(n 2 - p) 3 - An (n 2 - p) g + 8g 2 ) = . 

(3.15) 

Therefore, C 3 /A(27) is a complete intersection, given by a single (Calabi-Yau) hyper- 
surface in C 4 . 



3.2 Toric Varieties 

Having studied the first class of CY singularities, viz., the orbifolds, in some detail, 

let us move onto the next, and recently much- revived, class of geometries, the toric 

singularities. It turns out that here one can also write the Poincare series / explicitly, 

now in terms of the combinatorics of the given toric diagram D [37]. Mathematically, 

this is a nice extension of the Molien series. 

We first draw the graph dual 4 of the triangulation of D; this is the (p, g)-web [46], 

a skeleton of tri-valent vertices indexed by i G V. At each vertex i, the j-th (for 

incidentally, we remark that a convenient way of finding the dual (p, q)-web given a toric diagram 
is to take, for each pair of toric points, their cross-product, which then gives a vector perpendicular to 
the plane defined by the two said points. 
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j = 1,2,3) of the three coincident edges has charge with a a three-vector indexed 
by k, signifying the 3 charges. We remark that toric Calabi-Yau threefolds have three- 
dimensional toric diagrams whose endpoints are co-planar and this is why D and the dual 
(P) <?)-web are usually drawn on the plane. Here, we need to restore the full coordinates 
of the 3-dimensional toric diagram; whence, a has 3 components. With this notation, the 
Poincare series for D is (cf. Eq. (7.24-5) of [37] and also [47] for interesting mathematical 
perspectives): 

3 1 

P(h, t 2 , t 3 ; D) = J2U ^-^r • (3-16) 

Before we proceed, let us remark on the charges of coordinates t 1: t 2: t 3 . For toric 
varieties, we have three U(l) global symmetries: one is R-charge and the other two, 
flavor charges. In general, each coordinate tj is charged under all three U(l). For 
example, the R-charge of t 2 is given by the inner product of (0, 1, 0) and the Reeb 
Vector Vr = (bi,b 2 ,b 3 ). We recall that in the AdS/CFT correspondence the U(l) R- 
symmetry is dual to a special Killing vector, the so-called Reeb vector (cf. e.g. [32]), 



3 



which can be expanded as Vr = X^?ir> where 0j are the coordinates parametrising 



i=l 



'■d<j>i 



the T 3 -toric action. It is a very important quantity in toric geometry. 

It is possible to make coordinate transformation (ti,t 2 , £3) — > (ti,t 2 , t 3 ) such that each 
coordinate U is charged under one and only one U(l). However, the charges of these 
new coordinates U in general are not even rational numbers (for example the R-charge of 
dP 2 ), so it is not proper to use this new tj coordinate to do the Poincare series expansion, 
which must have integer powers. Furthermore, as we have seen in the example of C 3 , 
sometimes we want to find the generating function of only one U{\) charge, for example, 
the R-charge. To do so, we merely make the substitution (t 1? t 2 , t 3 ) — > (t a , t b , t c ) for given 
a, 6, c G Z> and the expression will be simplified considerably. In a lot of cases, the 
interesting U(l) is a linear combination of all three ?7(l)'s as we will see shortly. 

Returning to eq. (3.16), we have some immediate checks. First, we recall that all 
Abelian orbifolds of C 3 (including C 3 itself) are toric. For example, the toric diagram 
for C 3 is a triangle with vertices (1, 0, 0), (0, 1, 0) and (0, 0, 1). The dual graph, i.e., the 
(p, g)-web, has a single vertex, with three edges in the directions (1,0,0), (0,1,0) and 
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(0, 0, 1) respectively: 



Toric Diagram 



(1,0,0) 




& 

(0,0,1) (0,1,0) 



(0,0,1) 



(p,q)-web ' ' 



(1,0,0) 



Hence, 



P(t 1; t 2 ,t3;C 3 ) 



(l-*l)(l-*2)(l-*3) = ^ k ml ' 



(3.17) 



This is precisely the result eq. (2.2) obtained from conventional methods in §2. 

For a less trivial example, take the conifold C (cf. [37] as well as an earlier result in 
[36]). The toric diagram has 4 points, with coordinates 



A = (0,0,1), B 



;i,o,i), 



C 



;i,i,i), d = (0,1,1) 



(3.18) 



as shown in the center of Fig. 1. There are two triangulations, giving two (p, g)-webs 
upon dualising; the two are related by flop transitions. Of course, we need to prove the 
counting is independent of such choices. Indeed, as the conifold is the building block 
to all flops in toric varieties, if we show that / is the same for the two (p, g)-webs, this 
would be true for all toric diagrams and thus we would be at liberty to make any choice 
of (p, g)-web. First, take the left one, given by the two triangles ABD and BCD. This 
gives us 2 vertices, with (p, g)-charges of, respectively, {(0,1,0), (1,0,0), (—1,-1,1)} 
and {(0,-1,1), (-1,0,1), (1,1,-1)}. Thus eq. (3.16) gives us 

1 1 



P(x,y,z; C) 



X) 



;i-y) i 



yj v ' v ' \ x v 
The second trianglization is given by ACD and ABC, giving us 

1 1 

+ 



(3.19) 



§ (i-y) (i-l) (i-*) (i-S) 



(3.20) 



It is easy to see that the two expressions eq. (3.19) and eq. (3.20) are the same. Thus 
indeed the generating function is independent of how we triangulate and how the dual 
(p, g)-web is obtained [37]. 



17 



Trianglization 1 (°>U) (1,1,1) 
D C 



(0,-1,1) 



(1,0,0) 




(-1,0,1) 



(p,q)-web 1 



(0,1,0) 



Q- 



A~ B 
(0,0,1) (1,0,1) 

Toric Diagram 



Trianglization 2 



(0,-1,1) 



(1,0,0) 




(-1, 



(p,q)-web 2 



(0,1,0) 



Figure 1: The toric data for the conifold C. There are two triangulations, related to by 
flops, and thus two (p, g)-webs. We see in the text that they lead to the same counting. 

3.2.1 Refinement: £/(l)-charges and Multi-degrees 

We see from eq. (3.16) that for toric varieties the counting is more refined than the 
Molien series eq. (3.1) as the latter only counts invariants of total-degree. There seems 
to be a straight-forward generalisation. In order to get the number of single-trace GIO's 
given the R-charge of each field, the Molien counting seems to be refinable to counting 
the number of independent polynomials of a given multi-degree (ii, . . . ,i 3 ). This is done 
by generalising the Molien series to: 



M(t, G) 



\G\ det(I — diag(ti, 



Eh. . fi fk 



(3.21) 



The caveat is that now the coefficients &i lr ..,i fc are no longer guaranteed to be integers 
for general groups. This corresponds to the fact that the invariants are not monomial in 
general, but, rather, polynomial. For example, in eq. (3.6), at degree 6, there is a single 
invariant, built of the sum of two monomials, of multi-degree (5, 1) and (1,5), each of 
which is not an invariant. 

Therefore, this refinement only makes sense in case there is a corresponding con- 
served charge associated with a U(l) that is part of the isometry of the CY manifold. 
The isometry of C 3 is SU(4) with rank 3. There is thus a maximum of 3 charges cor- 
responding to the maximal subgroup of the isometry group of the CY manifold. If the 
manifold is toric then there is a T 3 fibration and therefore a total of three U(l) charges 
and the index would be a function of 3 variables. If the manifold is not toric then in 
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many cases the isometry group has a rank smaller than 3 and in most cases in fact is ab- 
sent. Nevertheless there is at least one charge, counting the R charge, that corresponds 
to the choice of complex structure of the manifold. 

To summarize, there are some cases in which the rank of the isometry group is 2 
and in most cases the rank is 1. All cases in which the rank is less than 3 are non-toric. 
An example for a manifold with rank 2 is the set of complete intersection manifolds 
of the form x 2 + y 2 + z 2 + w k = 0. Is has a clear 5*0(3) isometry acting on the first 
3 coordinates and together with the natural degree of the variables form the isometry 
group 577(2) xU(l). For the case k — 2 the isometry grows to 50(4) x U(l) and having 
rank 3 indeed confirms that the manifold is toric - the familiar conifold. We will study 
this geometry again in §5.3.2. 

An example of a manifold of rank 1 is any C 3 -orbifold with a full non-abelian sub- 
group T of SU(3). For T in SU(2), we still have Af = 2 SUSY, so the global isometry 
is 577(2) x U(l), of rank 2. Indeed though there is no refinement for 3 charges, there 
should still be a refinement of two charges since the rank is 2. The first charge will 
denote the Cartan charge of the 577(2). We could, for example, define a degree which 
counts how many rr's and y's together, and another degree which counts how many z's. 
This is the reason that in trying to implement the refinement on 77 4 we found factors 
of 1/2. There is simply no corresponding conserved charge which corresponds to this 
generalization. 

3.2.2 The Family 

An infinite family of toric CY 3-folds which has recently attracted much attention, 
because of the discovery of explicit CY metric thereon, is the Y p ' q, s (cf. e.g. ,[48, 49, 50]). 
Let us now do the counting for these. The toric data is given by O = (0,0, 1), A = 
(1,0, 1), B = (0,p, 1) and C = (— l,p — q, 1). We have drawn it at the left hand side 
of Fig. 2. As indicated, we take the trianglization as connecting the point T a = (0, a, 1) 
to A and C with a — 1, ...,p (so T p = 77). Now we have 2p triangles given by T a AT a+1 
and T a CT a+ i, a = 0, ...,p — 1. For triangle T a AT a+ i we have the following charges and 
corresponding term 



{(1,0,0), (a,l,-a), (_ a -1,-1, a + !)}=}► 



1 




(3.22) 
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B » (0.P.1) 



Ypq 



(-i.p-i.i) 




(O.p.l) 



Xpq 



(-l,p-q+l,l) 



(-l,p-q,l) 



(0,0,1) o 



A (7,0,iJ 




mo,7j O 



A (1,0,1) 



Figure 2: The toric diagrams for the spaces Y p ' q and X p,q . 
For triangle T a CT a+1 we have 

{(-1,0,0), (-a + (p - g), 1, -a), ({q - p) + a + 1, -1, a + 1)} 



(1-1) (l 



c l+a-p+<, z l + a 



(3.23) 



Putting these together we have 



p-i 

v,y,q; Y M ) = ]T 



^Z(l-x) (1-^) (l-^ 



1 — a ^l+a 



c l + a — p+g 2 : 1 + c 



(3.24) 

Knowing Y p ' q , it is easy to go on to X p,q . It differs therefrom by the addition of one 
point (—l,p — q + 1, 1) in the toric diagram. So we use the above trianglization of Y p,q , 
plus one more triangle given by (0,p, 1), (— l,p — q + 1, 1) and (— l,p — g, 1). This one 
gives the following vector and hence a new term to eq. (3.24): 

{(1,0,1), (q-l,-l,p),(-q,l,-p)} ' 



3.3 The Del Pezzo Family 



The last category of CY singularities widely studied in D-brane gauge theories is the 
cone over the 9 del Pezzo surfaces. These surfaces are P 2 blown up at n generic points; 
the cone is CY if n — 0, . . . , 8. There is a close cousin to this family, viz, the zeroth 
Hirzebruch surface Fq, which is simply P 1 x P 1 and the cone over which is also CY. 
It is well-known that for dP n=0t lj2)3 and for Fq, the space is actually toric (q.v. [11]). 
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Therefore, we can use eq. (3.16) to obtain the following: 



P(z,x,y;dP ) = 

P(z,x,y;F ) = 

P(z : x : y;dP 1 ) = 

P(z,x,y;dP 2 ) = 

P(z,x,y;dP 3 ) = 



+ 



d-x)(i-v) (i--) ^ (i_ S ) (1-^)^(1-1) (i_ f ) (1-4.) 



NKH^) + (i-iltHKi-f) ^ (i-x)(i-i)(i-^) ^ (i-i)(i-i)(i-^, 



+ 



+ 



+ 



(i— )(i-f)(i-f) " (1-^(1-^(1-^) " (i-*)(i-f)(i-^) "~ 



(1-1) (l-x)(l-f )(1-|) (1-1) (l-y)(l-^) 

1 1 1 1 



+ 



(1-,)(1-I)(1-^) - (!-*„*) 
1 1 1 



(l-i) (l-x»)(l-f) + (l-^) M (l-^) ^ 



+ 



^(i-x) (i-?/)(i-^) ^ (i--)(i-^)(i-^) + (i-i) (i-i) (i— y-) ' 
We include toric diagrams and the dual (p, g)-webs here for reference: 



(3.25) 




A=(-l -1,1) B=(1,0,1) 

c=(o,i,i) o=(o,o,i) 





(0,-1,0) 



(-1,2,1) 



(1,0,1) 






(6) dP 1 



(1,0,1) 




l (-1,0,0) (1,0,0) 



(0,-1,1) 

J 



1,1,-1) 



(-1,0,1) 



(0,1,0) i P4 pi 
(0,-1,0) J 



P3 



(IAD (-1,1,0)* 



(1,-1,0) 



V(-i,o,o))( 



"(-1,-1,1) 

(0,1,0) 



(0,1,1) (-1,1,1) 



(0,1,1) (-1,1,1) 



(c) dP 2 (d) dP 3 

For completeness, we also include the data for F as well as PdPA, the so-called 
pseudo dP 4 surface, first introduced in [51], which is obtained from blowing up a non- 
generic point of dP 3 so as to keep it a toric variety: 
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3.3.1 A General Formula for dP n 

We can see that taking the limit x = y — > 1 to relax the refinement in eq. (3.25) the 
expressions become very simple. In other words, we neglect the two U(l) charges carried 
by x, y and keep only the U(l) charge carried by t = z (note that this is not the R-charge 
but a linear combination of the three U(l) charges). The result counts the single-trace 
GIO's of a given total degree: 



f(t; dP ) = ±±^±|i, f(t; dPO = f(t; F ) = ^±^±|_, 
f(t;dP 2 ) = ^±f, f( t ;dP 3 ) = ^f 



(3.26) 



We shall see in the next section what it means to set x,yto\ and how all this relates to 
projecting 3-dimensional toric diagrams to 2-dimensions and to dimers. For the mean 
time, observing the pattern eq. (3.26) for the above 4 members of the del Pezzo family, 
we propose the following general expression for the generating function: 

f{t) {n) := /«,(*; dP n ) = 1 + ( ( 7 1 ~^ ) ) 3 + t2 , n = 0, . . . , 8 . (3.27) 

We remark that the result for F is the same as dP±. This is not surprising because they 
both, when having 1 more generic point blown-up, become dP 2 . Also, dPo is a Calabi- 
Yau over P 2 , it is in fact simply the orbifold C 3 /Z 3 which we will encounter again in 
§5.1.2. Furthermore, setting n = 4 gives agreement with the recent (P)dP^ [51] result 
of Eq. 5.29 of [27]. Indeed, we shall revisit the del Pezzo family, and give full credence 
to eq. (3.27) in §5.3.1. 
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4 Dimers, Toric Diagrams and Projections 



By now we have seen the Poincare series /oo in full action in counting single-trace GIO's. 
Before proceeding to finding the generating function for the multi-trace case, let us 
take a brief but important interlude in how the counting in / is pictorially realised for 
toric varieties. In due course, we shall see how the invariants emerge in slices of the 
3-dimensional toric cones and how such projections relate to dimers and 2-dimensional 
tilings. Indeed, it is these observations which initiated our original interest in this 
problem of counting GIO's. 

4.1 Example: Dimers and Lattices for C 3 

We begin by discussing the simplest toric CY 3-fold, C 3 , which was first mentioned in 
§2 and then in §3.2. Let us see how to represent the chiral ring in the dimer diagram 
of C 3 . The dimer for C 3 is well-known [16] and is drawn in Fig. 3. There is only one 
gauge group and it is represented by a hexagon. We recall the fundamental fact that in 
a dimer, the polygonal faces correspond to gauge groups, edges, (perpendicular) to fields 
and nodes, to superpotential terms. Thus, a BPS GIO in the chiral ring 5 can be thought 
of as a path from the origin to a polygon, and we shall show below that it is in fact path- 
independent. Since we consider only chiral operators, we represent the (holomorphic) 
operator by oriented lines crossing edges, such that when a line crosses an edge, the 
black vertex is to its left (recall that the coloring convention in a dimer has orientation 
built in). There are three holomorphic fields denoted by X, Y, Z. As mentioned before, 
the F-term relations here make these three operators mutually commutative. 

Now, we discuss the chiral GIO's in detail. We shall do so according to the number 
of levels. Here, we define level to mean the number of X, Y, Z fields inside the chiral 
operators. This was what we meant by degree in the aforementioned generating function. 
For clarity, we have enclosed each level with a dotted red circle in the diagram. At level 
1, there are only 3, given by Tr(X), Tr(Y) and Tr(Z). This has been shown in level one 
of Fig. 3. In the figure we have given also the 3 anti-chiral operators Tr(A), Tr(F) and 
Tr(Z) for reference. Henceforth, we shall use blue to denote the GIO's in which we are 
interested, viz., the chiral single-traces ones. The 3 here, of course, correspond to the 

5 Since we are studying BPS mesons, we may emphasize the relationship between the concept of the 
"extremal BPS meson" (introduced in [28]) and that of the "zig-zag path" in [21, 22, 24]. 
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Lattice of Chiral GIO at Level One Level Two Lattice Levcl Three LmicE 

Figure 3: The Dimer configurations and the lattice structure of GIO's for C 3 , exhibited at 
the first 3 levels. We have drawn some mixed chiral-antichiral GIO's as well for illustration, 
but the ones of our concern, viz., the chiral ones, are drawn in blue. 

3t term in eq. (3.2). 

Next, let us move to level 2 of Fig. 3. This time, we can cross two edges, as shown 
by the second red circle. A few remarks are at hand. First, for the hexagon denoted 
by the blue XY (which is inside the chiral ring), we have two paths to reach from the 
center. One is from the center to hexagon Y Z then to XY. Another one is from the 
center to XZ then to XY. The key point is that these two paths give the same element 
Ty(XY) in the chiral ring. So, our first conclusion is that for the chiral ring, a GIO 
depends only on the starting and ending point of the path in the dimer model and does 
not depend on the path itself. This in fact is generic for every dimer model and not just 
for the simple hexagonal model discussed here. See [21, 24, 28] for a proof of this point. 

Furthermore, we have in fact drawn not only the chiral ring, but also the anti-chiral 

ring and the mixed chiral-antichiral operators 6 at level 2. For the mixed operators, it 

6 We remark that the mixed operators are protected only in AT = 4 because of enhanced SUSY, in 
generic AT = 1 theories the "protected rings" are only the chiral or antichiral ones. The other 1/2-BPS 
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is easy to see that now the path in the dimer does matter. We start from the center, 
go up and then go down and get XX. Similarly we can go southwest to get YY and 
southeast to get ZZ. This means that we should put all three XX, YY and ZZ in the 
center hexagon. In another word, the one-to-one correspondence we found for chiral or 
anti-chiral ring is lost. Because of this complexity we will not discuss mixed operators 
further in this paper (the anti-chiral ring is isomorphic to the chiral ring and need not 
be addressed separately). 

Focusing on only the chiral ring we can see that there are 6 (chiral) GIO's at level 
2, corresponding to the term 6t 2 in eq. (3.2). Also, The result of level 3 is given at the 
right of Fig. 3. Here, we give only the chiral ring operators at proper hexagons in this 
figure. There is a total of 10 as shown, corresponding to the 10t 3 in eq. (3.2). 

4.1.1 Lattice Structure and Planar Slices 

We thus conclude that: 

A (chiral) GIO's at level n corresponds to a polygon in the dimer, which is a 
chiral-distance n away from the center. 

In the above, a chiral-distance is measured by segments of only chiral operators, i.e., 
black vertices to the left. We have drawn these chiral GIO's in blue in Fig. 3. With 
this one-to-one correspondence between chiral GIO's at a given level and hexagons in 
the dimer diagram, we can see that in fact we have a lattice structure in R 3 . Each 
integer lattice point (a, b, c) with a,b,c > corresponds to a chiral operator. The level 
of this operator is given by (a + b + c). In other words, level n is given by the plane 
perpendicular to vector (1,1,1) and has distance n from the origin. It is interesting to 
notice that the vector (1, 1, 1) is the Reeb vector of C 3 . Thus the degree we are counting 
is exactly the R-charge. We have drawn these plane slices for each level in our figure as 
well in Fig. 3. 

This lattice is something with which we are familiar! It is nothing other than the dual 
toric cone for C 3 . Indeed, the definition of a toric variety is that it is the affine spectrum 
of the ring of monomials obtained from raising the coordinates to the powers of the 
lattice generators, i.e., the invariant monomials. This is what we are doing above. Level 
1 gives the monomials which are obtained from the lattice generators of the cone; level 
protected operators like the currents do not form a ring. 
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2 gives the monomials obtained from the toric cone intersected with the (non-primitive) 
lattice points one further step away, etc. 

4.2 Example: Dimers and Lattices for the Conifold 

Next, let us discuss the conifold. Here, we will see explicitly how we must count the 
GIO's up to relations from F-terms, as was mentioned in the introduction. The toric 
diagram was given in Fig. 1 in §3.2. The dimer model is the brane diamond [52] drawn 
on T 2 [16] and is given in Fig. 4. There are two gauge groups so there will be two types 
of polygons which are labeled 1 and 2. It is easy to see that we can locate these 2 gauge 
groups at lattice points. More explicitly, if we draw the lifting of T 2 in M 2 , we can 
identify an integer lattice point (a, b) to gauge group 1 if a, b are integers or gauge group 
2 if a, b are half-integers. Now, since we are considering the single-trace mesonic GIO's, 
we can neglect gauge group 2 and consider the holomorphic paths connecting different 
lattice points of gauge group 1, i.e., integer lattice points in 2-dimensions. 

In this gauge theory, there are four bi-fundamental fields A 1 ,A 2 and Bi,B 2 . We 
define the following operators which are in the adjoint representation of gauge group 1: 

M ,i = AxSx, M 1>0 = A X B 2 , M_ lfi = A 2 B 1 , M„ _i = A 2 B 2 . (4.1) 

It is easy to check that the F-term relations tell us that all four commute and obey 
one non-trivial relation: 

M 0) iM _i = M lj0 M_ 1)0 . (4.2) 

We can map the above quantities into the dimer model. As we have shown above, 
the dimer model can be mapped to a 2-dimensional integer lattice. The operator Mo,i 
can be mapped to vector (0,1) so we can use it to connect points (0,0) and (0,1). 
Similarly, M 1)0 , M_ li0 , M _i map to vectors (1, 0), ( — 1, 0), (0, — 1), respectively. Using 
this mapping, a single-trace GIO is mapped to a path connecting point (0, 0) to (a, b) 
using the above four vectors. The non-trivial relation eq. (4.2) is nothing, but the 
statement that after following consecutively vectors (0,1) and (0,-1) (or (1,0) and 
(— 1, 0)) we come back to the starting point. 

Using this picture we can see the lattice structure of holomorphic GIO's. For level 
0, it is the origin (0, 0) and corresponds to the identity operator. For level one, we 
can use only one M itj to connect (0, 0) to nearby lattice points. Thus we have four of 
them Tr(M lj0 ), Tr(M_ 1;0 ), Tr(M ,i), and Tr(M _i). For level two, we need to use two 
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1 M(0,1) 


M(0,1)M(-1,0)» 


• M(0,1)M(1,0) 


M(-1,0)M(-1,0) 


M(1,0)M(1,0) 


• 




M(0,-1)M(-1,0) 


\M(0,-1)M(1,0) 


M(O-l) 


M(0,-1>. 




X 


M(i.0)M(-l.0) = M(0.1)M(0,-l) 



• M(0,1)M(0,1)M(1.0) 



M(0, 1 ) M( 1 ,0) M(- 1 ,0) = M(0, 1 ) M(0 - 1 )M(0, 1 ) 



Level Zero 



Level One 



Level Two 



Level Three 



Figure 4: (a) The dimer configuration for the conifold and the dual lattice structure of 
GIO's; (b) In more detail, the actual operators (with built-in relations) corresponding to the 
lattice points at the first 3 levels. 



Mi j operators. It is easy to get to lattice points (±2, 0), (±1, ±1) as well as (0, 0). For 
(0, 0) we have two ways Tr(M(0, 1)M(0, -1)) or Tr(M(l, 0)M(-1, 0)). But by relation 
eq. (4.2) they are the same so we should count only once. Similarly we can draw the 
level three lattice diagram as shown in Fig. 4. 

4.2.1 Planar Slices and Lattices 

Now let us find the 3-dimensional box which projects to the above 2-dimensional picture. 
The vectors (0, 0, 1), (1, 0, 1), (0, 1, 1) and (1, 1, 1) of the toric diagram of C generates 
an integral cone; we can find the generators of the dual cone to be v\ = (1,0,0), 
t>2 = (0, 1,0), t>3 = (0, —1, 1) and v± = (—1,0, 1). By definition, a lattice point in the 
dual cone is given by positive integer linear combinations of these four vectors. It is 
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special in our case that these four generators i>j have their endpoints co-planar 7 . It is 
easy to find the vector u orthogonal to the plane generated by v j as u = (1, 1,2). In 
fact, in this case, the vector u is precisely the Reeb vector, so the level we are counting 
is also the R-charge 8 . 

Now, we can see how this 3-dimensional lattice generated by Vi projects to the 2- 
dimensional lattice. For level one, it is given by four v iy since all of them have i>j • 
(1, 1, 2) = I. For the level two, we need to find vectors (x, y, z) such that (1) (x, y, z) = 
Ylt=i a i v i w hh a,i > and integer; (2) (x, y, z) ■ (1, 1, 2) = 2 which gives x + y + 2z = 2. 
From these conditions, we find the following 9 points: 

{(2, 0,0), (0,2,0), (0,-2,2), (-2, 0,2), (1,1,0), (1,-1,1), (0,0,1), (-1,-1,2), (-1,-1,2)} 

(4.3) 

which is exactly what we find in the dimer model. In general level n should have (n+ 1) 2 
points. 

Now let us check this using our Poincare series, which from eq. (3.19) is 

p(*,v,«c)=- (4,) 



To count the level, notice that the Reeb vector is u — (1, 1,2), which means that the 
R-charges of x, y, z are 1, 1, 2 respectively. In other words, we should replace x — > q,y — > 
q,z — > q 2 , yielding 

P(q; C ) = 7T^ = E(^ + 1 )V; (4.5) 

whereby giving us the required (n + l) 2 counting! 

In fact we can do better than that. Let us do the following replacement x — > xq, 
y — > yq and z — > q 2 . The expression is changed to 

xy(l — q 2 ) , 11, . 

- 1 + q(x + y + - + -) + ... . (4.6) 



(1 - qx)(l - qy)(q - x)(q - y) ' x y' 

Comparing this with toric data we can see that x, y represent the Cartan weight of 
SU(2) L x SU(2) R global symmetry for the conifold. More explicitly, for the two U(l) x 
U(l), x,y carry the charge of U(l) x = U(1) L + U(1) R and U(l) y = U(1) L - U(1) R . To 
see this let us consider, for example, M 0j i = AiBi. Because the (U(1)l, U(1) r ) charge of 

7 Note that though the toric diagram always has its vectors in a plane, as guaranteed by the CY 

condition, the dual cone is not so guaranteed. 

8 We remark that our Reeb vector differs in convention from that of [37]. Our (a, b, c) is their (c, a, b). 
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A 1 and B 1 is (|,0) and (0,|), we get immediately the (U(l) x ,U(l) y ) charge (1,0), i.e., 
the term x. Similarly, the terms y, |, ^ correspond to the operators M 0jl , M_ lj0 , M 0j _i, 
respectively. 

5 Counting Gauge Invariants: Plethystics, Multi- 
Trace and Syzygies 

We have now seen the generating function / which counts single-trace GIO's of a given 
choice of global charges for 3 large families of CY threefold singularities. What about 
the multi-trace GIO's? These are products of combinations of single-traces. We have 
called the generating function for counting these, g. We shall now see how g can be 
obtained from / using some nice combinatorics. We shall then see how the function 
which relates / and g has some remarkable geometrical properties as well. 



5.1 The Plethystic Exponential: From Single to Mult i- Trace 

Recall that in the above / should really be bacause we have taken the large N limit. 
Similarly, the quantity g we desire is really g^. Now, we showed in §2 that for C 3 , the 
relation between / and g is that of the plethystic exponential, PE (q.v. [53, 54]). 
This in fact holds in general: 

g(t) = PE[f{t)\ := exp j . ( 5 .1) 

Indeed, recalling eq. (2.6), we summarise the following relations, with the subscripts 
restored: 

9i = foe, /«,(*) = PE[h{t)l 9oo (t) = PE[ 9l (t)\ = PE[PE[h{t)\] . (5.2) 

We remark that, even for a list of variables ti=i,...,n, which are used in the refinement 
of counting discussed above, the expressions in eq. (5.1) and eq. (5.2) still hold, with 
obvious replacement. Namely, 



g(t u ...,t n ) = PE[f( tl , . . . , Q] := exp ^ /(^■■■,^ /(0,.--,0) ^ 



(5.3) 
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We can derive the statement eq. (5.1) explicitly by series-expansion. Let 

oo 
n=0 

be the Taylor expansion of the Poincare series = f(t). Thus, a n is the number of 
independent invariants at (total) degree n. Then, eq. (5.1) gives us 

(oo oo nk oo 1 \ 

n=0 k=l k=l / 

(1 \ 
- ^a„log(l-r)-ao^- . 
n=0 k=l J 

We see therefore that the /(0) term precisely regularises the sum and we obtain 

Pi^/(*)] =exp(-f>„log(l-r) ] =-55— ^ • ( 5 - 5 ) 

V n=l / Yl (! - tn ) a " 

n=l 

This expression is now in the standard Euler product form. Upon expansion of PE[f(t)], 
we would see that the coefficient for t m is the number of ways of partitioning m, each 
weighted by a n . This is precisely our required counting, i.e., the number of multi-trace 
GIO's at degree m. Hence, g(t) = PE[f(t)]. 

We have thus solved problems (1) and (2) posed in the introduction and have the 
generating functions / and g for large N. In fact, as before, we can refine our counting. 
In the above, we had a single variable t, a dummy variable associated with the total 
degree. Where permitted, as discussed in §3.2.1, we can have a set of variables U, one 
for each [/(l)-charge, and an associated multi-degree for these tuples of charges. In 
addition, we can introduce one more variable is, to be inserted into the summand. One 
could easily see that upon expansion, the power of v will actually count how many single- 
trace operators are present in each of the terms. In other words, for /oo(ii, • • ■ ,t m ) = 

oo 

E a P u-,pJT we have 

Pl,—,Pm=0 



(oo 
E 
k=l 

= ( n (i - ut T ■ • -C m ) api - 

\Pl,-,Pm 



/ o(tf,-,**>* 
k 



-1 



(5.6) 



note that due to the insertion of u, there is no longer a need to regulate the sum by the 
subtraction of /(0, . . . , 0). 
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5.1.1 The Plethystic Logarithm 

The inverse function of PE is also a fascinating one. It is called the plethystic logarithm 
[53]; one can in fact write it analytically: 

/(f) = PE-\g(t)) =jr^T MS^)) , (5-7) 

k=l 

where /i(k) is the Mobius function 

k has one or more repeated prime factors 

//(//) = { 1 k = 1 • (5.8) 

(— l) n A; is a product of n distinct primes 

As = PE[gi], so too does one have the relation = PE[fi]. Since our basic 
generating function is the Poincare series f = foo, for which we have had explicit results 
in §3, it is more convenient to write 

h = PE-'Uoo) . (5.9) 

One may ask what this function fi, which we briefly encountered in eq. (2.5), signifies. 
It has a remarkable geometrical property! 

The plethystic logarithm of the Poincare series, is a generating series for the 
relations and syzygies of the variety! 

We exemplify this statement with our familiar example of the Valentiner group from 
§3.1.2 in the next subsection. 

5.1.2 Plethystic Logarithm and Syzygies 

Using eq. (5.7), and recalling the Poincare (Molien) series / for A(27) from eq. (3.13), 
we see that 

-1 +f 3 -f 6> 

:-i+* 3 ) s 

The RHS terminates and is a polynomial! It is to be interpreted thus: there are 2 degree 
3 invariants, 1 degree 6 and 1 degree 9 invariant, these 4 invariants obey a single relation 
of total degree 18. Upon inspecting eq. (3.14) and eq. (3.15), we see that this is indeed 
the definition of C 3 /A(27) as a variety! 



fl = PE- 1 ( i-Ti _ I ) = 2t 3 + t 6 + t 9 - t 18 . (5.10) 
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Uit) = M(t; Z 3 ) = -± ^3- , (5.11) 



Now, C 3 /A(27), as a hypersurface in C 4 , is a complete intersection affine variety (i.e., 
the number of equations is equal to the codimension of the variety in the embedding 
space). How does the above work for non-complete intersections? We have an example 
readily available: the famous C 3 /Z 3 = (9 P 2(— 3) orbifold. In fact, being an Abelian 
orbifold, this is also toric and furthermore, it is also dP , being a cone over P 2 . So we 
have 3 ways to compute its Molien series from §3. Let us use the Molien series. The 
action is (x, y, z) — > ojz{x, y, z) and we immediately get 

1 + 7 1 3 + t 6 

(1-* 3 ) E 
whereby 

f\{t) = PE-^f^t)] = 10/ 3 -27t 6 +105t 9 -540/ 12 +3024£ 15 -17325/ 18 +C>(t 21 ) . (5.12) 

This is again in accordance with known facts! The equation for this orbifold is 27 
quadrics in C 10 , i.e., 10 degree 3 invariants satisfying 27 relations of degree 2x3 = 6 
(q.v. [38, 43]). We can determine these as follows. The 10 invariants are 

f 3 2 232 2 2 231 /r i o\ 

yi,...,w = {x , x y, xy , y , x z, xyz, y z, xz , yz , z } , (5.13) 

obeying the 27 quadrics 

{2/2 2 - 2/12/3, 1/22/3 - 2/12/4, 2/3 2 - 2/22/4, 2/22/5 - 2/12/6, 2/32/5 - 2/12/7, 2/22/6 - 2/12/7, 
2/42/5 - 2/22/7, 2/32/6 - 2/22/7, 2/42/6 - 2/32/7, 2/5 2 - 2/12/8, 2/52/6 - 2/12/9, 2/22/8 - 2/12/9, 
2/6 2 - 2/22/9, 2/52/7 - 2/22/9, 2/32/8 - 2/22/9, 2/62/7 - 2/32/9, 2/42/8 - 2/32/9, 2/7 2 - 2/42/9, 
2/52/8 - 2/12/10, 2/62/8 - 2/22/10, 2/52/9 - 2/22/10, 2/72/8 - 2/32/10, 2/62/9 - 2/32/10, 2/72/9 - 2/42/10, 
2/8 2 - 2/52/10, 2/82/9 - 2/62/10, 2/9 2 - 2/72/io} • 

(5.14) 

Therefore, eq. (5.13) are the 10 primitive invariants of degree 3, obeying 27 syzygies of 
degree 6 given by eq. (5.14). According to our rule, this should read 10t 3 — 27t 6 . These 
are precisely the first two terms of eq. (5.12)! Indeed, because we no longer have a 
complete intersection, the plethystic logarithm of the Poincare series is not polynomial 
and continues ad infinitum. What about the 105 and higher terms then, do they mean 
anything? We will explain this in §5.2. 

As a final example of the more subtle case of non-complete-intersection varieties, let 
us take the C 3 /Z 5 orbifold, with action (x,y,z) — > (u 5 x, cufy, w\z). We obtain: 
— 1 -4- 1 — 3 / 3 -I- / 4 — 3 / 5 -I- t 7 — t 8 

M(t;Z 5 ) = + f L = l +3t 3 + 2t 4 7t 5 + 5t 6 + 4t 7 1U 8 

(-l+/) 3 (l+/ + / 2 +t 3 + t 4 ) 2 

+9t 9 + 18t 10 + 15 1 11 + 13t 12 + 24 1 13 + 21 1 14 + 34 1 15 + C(t) 16 , (5.15) 
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giving us 

PE'^f^t)] = 3t 3 +2t 4 +7t 5 -t 6 -2t 7 -13t 8 -12t 9 -14t 10 +14t 11 +34t 12 +72t 13 +47t 14 +C(t) 

(5.16) 

We can find that the 3 invariants of degree 3, 2 of degree 4, 7 of degree 5 are 

2/i 12 : = {xy 2 ,xy z,xz 2 ,x 3 y,x 3 z,x 5 ,y 5 ,y 4 z,y 3 z 2 ,y 2 z 3 ,y z*,z 5 } . (5.17) 

We can easily find, using Grobner algorithms, all relations amongst these 12 invariants, 
giving us 1 in degree 6 (a quadric in the 3 degree 3 invariants), 2 in degree 7, 13 in 
degree 8, 12 in deree 9 and 16 in degree 10. All this is in almost in exact agreement 
with eq. (5.16), with the only exception being that there are 16 degree 10 relations and 
not 14. Together with the issue of the higher terms in the C 3 /Z 3 case, we now address 
this discrepancy in the next subsection. 

5.2 Plethystics: A Synthetic Approach 

We have now witnessed the astounding power of plethystics in the counting problem 
and have moreover noted a tantalising fact about the geometry of the variety and the 
(plethystic logarithm of) the generating function for the GIO's in the gauge theory. Let 
us now attempt to argue why some of the above examples should work. First we note 
that the Poincare series /, when finally collected and simplified, is always a rational 
function. In particular it has a denominator of the form of products of (1 — t k ) with 
possible repeats for k; the numerator is some complicated polynomial. We will call this 
the Euler form. The point is that the coefficient in front of the t k is always unity and 
we conjecture that this is a property of the Poincare series of concern 9 . 

When we are taking the plethystic logarithm of /, due to the explicit expression 
eq. (5.5), we are trying to solve the following algebraic problem: find integers b n such 
that 

/(*) = ^o-^ , (5-18) 

n (i - t n ) b - 

71=1 

OO 

where f(t) is a given rational function in Euler form. Note that P-E 1-1 [/(£)] = b n t n , 

n=l 

unlike the Poincare series herself, need not have all positive b n . Because f(t) has Euler 

9 In fact, all Poincare series we have encountered, orbifold, toric, etc., have this property. We do not 
have a rigorous proof of this right now and leave it to the mathematically inclined. 
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form, the denominator of eq. (5.18) is immediately taken care of. In other words, because 
f(t) has denominator in the form of products of (1 — t k ), all positive values of n and b n 
are just read off. These are low values of n and correspond, in the Molien case, to some 
of the small invariants, including the primitive ones. However, there is still a numerator 
in f(t), often of complicated form. This will give negative b n contributions to the RHS 
of eq. (5.18), which correspond to the relations. 

Take A(27) as an example. We need to find b n such that 

1 - 1 3 + 1 6 (1-t 18 ) 1 



a-* 3 ) 3 "(i-t 6 )(i-t 9 )(i-t 3 ) 2 n(i-^) 6 " 

n=l 

where we have used the identity 

(i - 1 3 ) a - 1 18 ) 



(5.19) 



1 - t 3 + t 6 . (5.20) 



This rational identity is crucial and expresses even the numerator of / into Euler form. 
Now we can read out the solution: the denominator contributes terms +2t 3 , +t 6 and 
+t 9 while the numerator contributes -t 18 . Thus PE- l [M{t)} = 2t 3 + t % + t 9 - t 18 . In 
other words, there should be 2 degree 3 invariants, 1 each of degrees 6 and 9, obeying a 
single relation of degree 18. The fact that the numerator can be factorised into (finite 
polynomial) Euler form dictates that the plethystic logarithm of / terminates in series 
expansion as was seen in eq. (5.10). The moral of the story is that 

Findings relations in this language corresponds to finding appropriate factorisa- 
tions of the numerator into Euler form. 

Of course, not all Poincare series have polynomial plethystic logarithms. This is just 
the statement that not all polynomials afford identities of the type eq. (5.20). In general, 
the product on the RHS of eq. (5.18) must be infinite to accommodate those which 
cannot be put into finite Euler form. These correspond to non-complete intersection 
varieties. Take C 3 /Z 3 , we have 

l + 7t 3 + t 6 

Indeed, no rational identity can express 1 + 7 1 3 + 1 6 in finite Euler form so the plethystic 
logarithm will not terminate. Now, as promised earlier, we can explain the higher terms 
such as the 105 and 540. In this example, there are 10 basic invariants and there are 27 



M(t;C 3 )= Z Zs ■ ( 5 - 21 ) 
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relations amongst them. This explains the first 2 terms in eq. (5.12). This is seen above 
because if we were to write 1 + 7t 3 + 1 6 in Euler form, we would obtain 

(\ _ f6\27/i _ f 12\540 

thus we get the +10t 3 term from the (1 — t 3 ) 3+r factor in the denominator and the — 27t 6 
term from the (1 — t 6 ) 27 factor in the numerator. 

Of course, one finds the 27 relations by finding syzygies among the 10 primitive 
invariants. The reason we can do this is of course Theorem 3.1 which dictates that we 
need not go beyond degree \G\ = 3 to find all basic invariants which generate the entire 
invariant polynomial ring. Instead, if we found the syzygies for the entire invariant ring, 
we would get the higher terms. That is, we should, considering the expansion of the 
Molien series 

1 + 7t + * =l + 10t 3 + 28t 6 + 55t 9 + 91t 12 + 136t 15 + 190t 18 + ... , (5.23) 
(1 - t 3 ) 3 

consider all 10 + 28 + 55 + . . . invariants as polynomials in 3 variables and find all 
their syzygies; this should give the higher terms. Of course this cannot be done all at 
once, but nevertheless we can consider the process stepwise: first, syzygies for 10 of 
them, then 10 + 28 of them, etc. In principle, if we only wish to know up to some 
degree, we only need to find syzygies for invariants up to that degree. This algorithm 
is the precise analogue of the infinite product expansion of eq. (5.22) into Euler form, 
which serves as successive approximation to the LHS in eq. (5.22). This also explains the 
discrepancies in the case of C 3 /Z 5 as seen above. In these cases where the Euler product 
is non-terminating, and rational identities become infinite products, the syzygies should 
thereby receive stepwise corrections. We should be able to arrive at the right answer 
after some finite number of steps if we only wish to know the terms up to a desired 
order. 

Let us check up to second order in this example of C 3 /Z 3 by finding the syzygies 
amongst the 10 basic invariants of degree 3 and 28 degree 6 invariants. We find, using 
[42], 595 relations: 55 of degree 6, 225 of degree 9 and 315 of degree 12. This thus reads 

10t 3 + 28t 6 - 55t 6 - 225t 9 - 315t 12 = 10t 3 - 27t 6 - 225t 9 - 315t 12 . (5.24) 

Good, we reproduce the first 2 terms of Pi?~ 1 [M(t)] and have the next 2 terms. This is 
only up to order 2, i.e., finding syzygies among 38 polynomials. At next order, we would 
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have to find relations among 10 + 28 + 55 = 93 polynomials and correct the t 9 and t 12 
coefficients; the computation becomes increasingly strenuous for the computer 10 . 

5.3 Complete Intersections 

We see from the preceeding arguments that the most powerful avatar of the intimate re- 
lations between plethystics and syzygies is realised in complete intersections, especially 
in single hypersurfaces. We have seen that A(27) is one such example of the hypersur- 
face. The key feature for this class of varieties is that the series for f\ = P£'~ 1 [/ 00 (t)] 
terminates and is polynomial. This is nice because if we knew the defining equations 
and the degrees of the various pieces, we could re-construct the Poincare series and find 
the number of invariants in the gauge theory! This is independent of whether the variety 
is orb if old or toric, but should hold in general. In fact, we do not even need to know 
what the gauge theory is! We shall see, in §6.2, an inverse application of this philosophy, 
where we shall construct a variety with desired gauge invariants. 

5.3.1 del Pezzo Family Revisited 

Take a non-orbifold, non-toric, single hypersurface, the famous cubic in P 3 ; this is the 
cone over the 6-th del Pezzo surface. From eq. (3.27) and eq. (5.7), we have 

/(*; dP 6 ) = l * t _ + J =► PE-'Uit; dPs)} = At - t 3 , (5.25) 
(1 t) 

which says that there should be 4 linear invariants, obeying 1 cubic relation; precisely 
the definition of dP 6 . 

Another illustrative example is dP$; here we shall go beyond projective spaces, but 
rather to weighted projective spaces. We shall see that the plethystic logarithm still 
works. We know (cf. [55]) that dP 8 as a surface is given by a single equation in WPf 12 3 - 
Again, from eq. (3.27) and eq. (5.7), we find 

/(*; dP 8 ) = PE'Vit; dP 8 )\ = 2 1 + t 2 + t 3 - t & . (5.26) 

(1 t) 

This is again correct: 2 degree 1, 1 degree 2 and 1 degree 3, obeying a single degree 6 

relation. This can only happen in a weighted projective space, viz., W¥\ 12 ^- Thus, our 

10 This alternative addition of invariants and substraction of relations is reminiscent of the characters 
of minimal models and the removal of null-states using the Kac determinant in 2-dimensional conformal 
field theories. 
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proposal eq. (3.27) is again confirmed. We note that, upon comparing eq. (5.26) and 
the result eq. (5.10) for A(27), the /-functions are the same, only with the replacement 
t — > t 3 . This does not surprise us, indeed (cf. e.g. [7, 56, 57]) A (27) is known to be a 
special point in the moduli space of c?P 8 's. 

In fact, all del Pezzo surfaces for n > 4 are complete intersections (cf. e.g., eq.(3.2) 
of [55] and also [58]). We check against eq. (3.27) and find complete agreement. For 
clarity, let us tabulate these results: 



dP n 


/ = /oo(f) 


fi = PE-Wt)] 


Defining Equations 


5 


l+2t+t 2 
(I-*) 3 


ht-2t 2 


2 degree 2 equations in P 4 


6 


l+t+t 2 
(l-*) 3 


At-t 3 


1 degree 3 equation in P 3 


7 


1+t 2 
(I-*) 3 


3t + t 2 -t 4 


1 degree 4 equation in WTP 3 lil;2 


8 


1-t+t 2 
(I-*) 3 


2t + t 2 + t 3 -t 6 


1 degree 6 equation in WFf 1 2 3 



Therefore, for the entire del Pezzo family, members to 3 are checked by toric meth- 
ods while 5 to 8 are complete intersections. The only one remaining is dP 4 and from 
eq. (3.27), 

X _i_ 3 1 _i_ f 2 

f(t;dP A ) = — = l+6t+16t 2 +3H 3 +5H 4 +76t 5 +106t 6 +14H 7 +181t 8 +226t 9 +C(t 10 ) 

(1 — t) 

(5.28) 

predicts the single-trace GIO counting for this variety. The equation for this variety 
[51, 58] is the (non-complete) intersection of 5 quardrics in P 5 (cf. also eq. 5.29 of [27]). 
Expanding the plethystic logarithm of / in this case gives 

PE^[f(t; dP 4 )} = 6 1 - 5 1 2 + 5 1 3 - 10 1 4 + 24 1 5 - 55 1 6 + 120 1 7 - 0{f) (5.29) 

We see that the first 2 terms are actually correct: there are 5 degree 2 relations in 6 
variables! 

For full reference, we tabulate below the other members of the del Pezzo family, 
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these are non-complete intersections: 



dP n 


/ = /oo(f) 


h=PE~ l [f(t)} 


Defining Equations 





l+7t+t 2 
(1-tf 


10 1 - 27 t 2 + 105 i 3 - 540 1 4 + 3024 1 5 + C(t) 6 


(10|2 27 ) 


1 


l+6t+t 2 
(1-tf 


9 1 - 20 1 2 + 64 1 3 - 280 1 4 + 1344 1 5 + 0{tf 


(9|2 20 ) 


2 


l+5t+t 2 
(1-tf 


8 t - 14 t 2 + 35 1 3 - 126 1 4 + 504 1 5 + 0(tf 


(8|2 14 ) 


3 


1+U+t 2 
(1-tf 


7t - 9 1 2 + 16 t 3 - 45 1 4 + 144 1 5 + 0{tf 


(7|2») 


4 


l+3t+t 2 
(l-tf 


Qt-5t 2 + 5t 3 - 10t 4 + 24t 5 + 0{tf 


(6|2 5 ) 



(5.30) 



Here, we have computed these defining equations using fat-point methods on P 2 [42, 59]. 
We have used the notation, in the above table, that (m\pf . . .p q k k ) denotes q 1 equations 
of degree p±, q2 equations of degree P2, etc., all in m variables. The first memeber, dPo, 
is of course C 3 /Z 3 as studied in detail in eq. (5.13) and eq. (5.14). Furthermore, as 
mentioned when we first presented eq. (3.27), F has the same f(t) as dP\. Indeed, we 
can study the degree 2 Veronese-Segre embedding of P 1 x P 1 into P 8 and see that F also 
has defining equation in 9 variables as (9|2 20 ). The precise forms of these equations, of 
course, differ from those of dP\. Thus, dP\ and F are in different points of a complex 
structure moduli space. 

To compare and contrast, we include the f\ results for the ADE-series addressed in 
§3.1.1; indeed these are all complete intersections - in fact, single hypersurfaces - so /i, 
the plethystic logarithm of the the Molien series should be polynomial: 



G C SU(2) 


f 1 = PE~ 1 [M(t; G)] 


Defining Equation in C[u, v,w] 




t 2 + 2t n - t 2n 


uv = w n 


Dn+2 


f 4 + f 2n + f 2n+2 _ f 4n+4 


u 2 + v 2 w = w n+1 




f6 + t 8 + t 12 _ t 24 


u 2 + v 3 + w A = 


E 7 


f 8 + f 12 + f 18 _ t 36 


u 2 + v 3 + v w 3 = 


E 8 


t 12 + t 20 + t so _ f m 


u 2 + v 3 + w 5 = 



(5.31) 



5.3.2 Example: The Hypersurface x 2 + y 2 + z 2 + w k = 

Now, let us try another family of complete intersection 3-folds, viz., x 2 + y 2 + z 2 + w k = 
in C 4 . For k — 1, this is just C 3 , for k — 2, it is the conifold C. For k > 2, the theory 
is studied in [14]. However, for k > 3, [60, 61] recently showed that there is no Sasaki- 
Einstein metric, whereby making the AdS/CFT correspondence a little ambiguous here. 
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For k = 2n even, we have x, y, z being degree n and w being degree 1. From this we 
can read out fi — t + 3t n — t 2n . Thus we can calculate that 



Zoo (A; = 2n) = PE[t + 3t n -t 



2„1_ (I-* 2 ") 



(1 -t)(l -t") 3 ' 

Similarly, for /c odd, we have x,y,z being degree A; and u>, degree 2. From this we can 
read out that f\=t 2 + 3t k — t 2k and whence 

foo(k) = PE[t 2 + 3t k - t 2k ] = (1 _^ 2 ~f^ tfc) 3 , , k odd . 

To demonstrate, we list the series expansion of the Poincare series for k = 1 to k = 5 
and we find 

/oo(l) = /«>(*; C 3 ) = 1 + 3t + 6t 2 + 10t 3 + 15t 4 + 21t 5 + 28t 6 + 36t 7 + ... 

/oo(2) = /«>(*; C) = 1 + At + 9t 2 + 16t 3 + 25t 4 + 36t 5 + 49t 6 + 64t 7 + ... 

7^(3) = l + t 2 + 3t 3 + t 4 + 3t 5 + 6t 6 + 3t 7 + 6t 8 + 10t 9 + 6t 10 + ... 

/oo(4) = l + t + 4t 2 + 4t 3 + 9t 4 + 9t 5 + 16t 6 + 16t 7 + 25t 8 + 25t 9 + 36t 10 + 36t n + ... 

7^(5) = l + t 2 + t* + 3t 5 + t 6 + 3t 7 + t 8 + 3t 9 + 6t w + 3t u + 6t 12 + 3t 13 + 6t u + ... 

5.4 Refined Relations: The Conifold Revised 

In all of the above, we have used the generating function with a single variable t. How 
does this all work if everything is refined fully so as to contain a tuple of dummy variables 
for the various [/(l)-charges? We shall now see that the relations are still explicitly 
encoded by f x . We mentioned early on that the F-term relations are automatically 
built into the counting. Indeed, for C 3 , we do not have these relations - just that 
x, y, z commute. For the conifold, we have the simplest demonstration that fi contains 
relations. 

Recall the expression for the Poincare series in eq. (4.6). Now, let us take the 
multi-variate plethystic logarithm [53], which for f 00 (t 1 , . . . ,t m ) is, recollecting eq. (5.6) 

PE-'iUt,, t m )\ = log(/oo(ff, • • • , t k J) . (5.32) 

k=i 

The result is 

^ PE ^(1 -qx){l - qy)(q - x)(q - y)^ x + 9% + y + 9V 9 
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Indeed, fi is polynomial because the conifold is complete intersection. There are four 
invariants, corresponding to qx, -, qy, -. If we would write just these generators without 

x y 

any subtractions, this would be merely the result for C 4 . Therefore it is not enough. 
To put the relation, we notice that (qx)(q/x) = (qy)(q/y) = q 2 and therefore we should 
subtract one combination of q 2 . We thus reproduce eq. (5.33). The procedure is simple 
and analogous for complete intersections. However, for non-complete intersections, once 
we make the subtraction we are taking away too much. We must therefore compensate 
by adding those which are subtracted, etc. ad infinitum, just like the non-terminating 
series explained in §5.2. 



6 Asymptotics and the Meinardus Theorem 

We have encountered, in the preceeding discussions, many infinite products of Euler 
type. Indeed, we recall that the generating function for multi-trace GIO's is 

oo oo 

g(t) := = PE[f(t)\ = for f(t) = £ a m t m . (6.1) 

n=0 Y\ (1 - t n Y n m=0 



n=l 



In the case that all a n = 1, g(t) is the Euler function, or, up to a factor of t~24, the 
Dedekind 77-function. This is our familiar generating function for the number of ways of 
partitioning integers. The Hardy- Ramanuj an equation gives the asymptotic behaviour 
of p n and was what gave rise to the Hagedorn temperature (q.v. [64]). It is, needless to 
say, important to find analogous asymptotic behaviours for general a n . This would give 
micro-state counting for our quiver gauge theories. 

Luckily, this generalisation of Hardy- Ramanuj an is known. This is a result due to 
G. Meinardus [62] (q.v. [63], to whose notation we adhere, for some explicit results). 
Meinardus' theorem states that the asymptotic behaviour of p n in eq. (6.1) is: 



1 

n s+r(i + _) (AT (a + 1)C(« + 1)) ; 
a 



[l + 0(n- C3 )), (6.2) 



p n ~ Cin C2 exp 

if the Dirichlet series for the coefficients a m of /, defined as 

00 

D(s) := % Ms) > a > , (6.3) 



rrr 

m=l 



converges and is analytically continuable into the strip —Co < Re(s) < a for some 
real constant < C < 1 and such that in this strip, D(s) has only 1 simple pole at 
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s = a G M + with residue A. The constants in eq. (6.2) are 

d = e fl '(°> / Ml> + l)C(a + l))^gf , 
a/27t(q; + 1) 

D(0) - 1 - f 

C 2 = (6.4) 
and C 3 some positive constant. 

6.1 Example: C and Dedekind rj 

For example, when all a m — 1, we have the usual partition of integers and the Dirichlet 

00 

series is just the Riemann (^-function. The generating function f(t)— ^2 t m is of course 

simply and the geometry is that of the complex line C. Using the above results of 
Meinardus, we have 

a=l, A=l, D(0) = ~, D'(0)=e-^^; C x = -±=, C 2 = -1 , 
giving us 

Pn ~ ^^(1 + 0(n" c )) , 
precisely the Hardy- Ramanuj an result. 

6.2 Example: MacMahon Function and a Riemann Surface 

Next, consider 

P , v 1-t + t 2 1 

u tj n(i-r)" 

71=1 

As P-E -1 [/(£)] = t + 1 2 + t 3 — t 6 , this is a complete intersection, given as a hypersurface 
of degree 6 in W¥l 23 . The dimension is therefore 1 and is hence a Riemann surface. 
The total space is an affine cone over this surface and is of dimension 2. We can embed 
WFI23, through a Veronese- Segre map, into P 6 and see that the genus of the Riemann 
surface is 1. Alternatively, we can projectivise the weight-one coordinate in IUP123 and 
simply obtain a standard elliptic curve. Therefore, the geometry is a cone over a torus! 
The generating function g(t) is the well-known MacMahon function [65]. We see that 

a = 2, A=l, D(Q) = --, D'(0) = --log(G,); d = > C 2 = --, 

1Z LZ 236 G; V-J7T OO 
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where Gi := ^ ~ 1.28243 is the Glaisher constant. Hence, 



p n ~ -it ^ y n 36 exp -^n3 ) . (6.5) 
2se Gi V3 7T y 2s 

6.3 Example: Our Familiar C 3 

Returning to something we have encountered earlier, let us attack the C 3 example of 
eq. (2.11). We can now find the coefficients therein! Using eq. (5.5) we have that 

g{t; C 3 ) = PE[f(t; C 3 )] = , a n = \{n + \){n + 2) . (6.6) 

n(i -#»)«« 

n=l 

We can readily see that D(s) = \ (C(-2 + s) + 3C(-1 + s) + 2((s)). We see that there 
are 3 poles, at 1,2 and 3. Of course, Meinardus Theorem requires that there be only one 
pole within a strip. Thus, one must consider one monomial of a n at a time and consider 
the break-down 

9(t; C 3 ) = := gi (t)g 2 (t)g 3 (t), 

n (i - t n ) an i n (i - tn Y n2 n (! - tn ) ans 

n=l n=l n=l 

with a ni = \r? , a n , 2 = |n, a„ 3 = 1. Applying Meinardus and defining <?j=i,2,3 : = 

oo 

^2 Pi(n)t n , we have that 

n=0 

C'(-2) 

-n s exp ( r Ji4 I , Jh\n) ~ — r 



3 2 



/ \ e LL ^r 1 _S ^22**- 3\ , s eB((3)^ _1Z /3(3C(3))3 

Pi{n) ~ -^-5 — r« 8 exp p 2 ( n ) ~ -^^-—n 24 exp ^ ^ ;; 

225 155 \3153 / G^v^ V 

Ps(n) ~ —^r n ~^ exp ^-^n^ j . 

(6.7) 

Therefore, we have the convolution p(n) = Piij-)p2{s)p^{t) and since the expo- 

r+s+t=n 

nential growth of p\{n) dominates over the other two, for large n 

p{n) ~ Pi(n) . (6.8) 



7 Single- Trace and Multi- Trace for Finite N 

We have, in all preceeding discussions, made the important simplification of taking N, 
the matrix size of the operators, to infinity, whereby decoupling spurious relations which 
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arise from the lack of commutativity among the matrices of finite size. As mentioned 
in the introduction, the problem of counting for finite N is a significantly more difficult 
one. Nevertheless, we shall see in this section that the plethystics are still applicable. 

We consider the problem of counting BPS states of Af = 1 supersymmetric quiver 
gauge theories for N finite; N is the number of D3-branes at the tip of the CY cone. 
We denote the generating function for multi-trace GIO's by g N . This problem is of 
significant interest, for instance, to studying phase transitions and AdS$ black holes. 
We already considered in the previous sections the functions g\ = /oo and = PE[gi], 
and we are going to propose that it is still quite simple to reconstruct g^ in terms of g\. 

oo 

Suppose the single-trace generating function is given by gi(t) = fooit) = a nt n , 

n=0 

then we can construct the following function: 



g^-t):=T\- (7.1) 

m=0 y ' 

We immediately notice a strong similarity to the Euler product form of the plethystic 
exponential introduced in eq. (5.1), eq. (5.5) and especially eq. (5.6). 

We now propose that the finite N multi-trace generating function g^if) is simply 
given by the expansion 

oo 

Y,9N(t)u N = g{u;t) . (7.2) 

N=0 

We have 2 limiting cases to check, viz., g\ and g^, with which we are now quite familiar. 
First, we note that 

igfrt)^ k n ( i- y ^ = ^ f )E(rr^)- 

fc=0 v ' m=0;m^k v ' k=0 v ' 

Furthermore, since g{0;t) = 1, we have that d u g(0;t) = YlT=o a k^ k ■ Therefore, the 
coefficient of v in eq. (7.2) is indeed our g\. 

d„g(v,q)\v=o = 9i(q) ■ (7.3) 

Next, let us check whether the iV-th coefficient for N — > oo gives our g^. This 
coefficient can be found by considering the limit 11 lim(l — v) aa g(y\ t) which extracts the 

V— »1 

11 In this paper we always have do = 1, since the only operator of vanishing scaling dimension is the 
identity. We can see this explicitly in all the examples we have given throughout. 
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large iV-term in the series expansion. We see that 



oo 



1 



hm(l-u) ao g(u;t)= T] 



(1 - t m ) 



PE[ 9l (t)} . 



(7.4) 



m=l 



Therefore, our expansion eq. (7.2) has the property that its large N coefficient is the 
PE of the linear coefficient, precisely what is required of g^. We will see in the ensuing 
text why the expansion does what it is supposed to. 

It is very interesting to compare eq. (5.6) and eq. (7.1). From it we can see that the 
parameter v in eq. (7.1) can have two different interpretations: 

I. It counts the number of single-trace GIO's in a multi-trace GIO for the limit of 
matrix rank N — > oo (here we include the single trace of identity as well) as in 



II. It counts the number of multi-trace GIO's for matrix rank N given by the finite 
power of v as in eq. (7.2). 

Naively these two counting problems seem to be unrelated, but our proposed formula 
eq. (7.1) indicates that they are the same. 

Formulae (7.1) and (7.2) give the general solution for counting multi-trace BPS 
GIO's, for a finite number N of D3-branes. In fact, the relation between and gw, in 
general, still obeys the plethytic exponential as was in eq. (2.8) and eq. (5.1), which we 
summarise now (for a list of variables ij): 



Symmetric Products and Moduli Spaces: We can in fact re-examine the finite 
N counting from another perspective. The standard lore for N D3-branes probing a CY 
manifold, X, is that the moduli space of vacua, Ai vac (N; X), is the symmetric product 
of iV copies of the CY manifold, 



where Sn is the permutation group of N elements. Following our general line in this 
paper we can now state two important relations: 



eq. (5.6); 




(7.5) 



N 



M vac (N;X)=S N (X): 



(7.6) 
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1. Qn counts multi trace operators for one D3-brane on the symmetric product of iV 
CY manifolds A4 vac (JV; X): 

X N X N 
g N (t; X) = 9l (t; — ) = — ). (7.7) 

on djv 

Alternatively we can think of it as the Poincare series for M vac (N; X). Further- 
more, from the second equality we conclude that gjy also counts the single trace 
operators on Ai vajC (N; X) in the limit in which there are no matrix relations at 
all, N -> oo. 

2. /jv counts single trace operators for one D3-brane on the symmetric product of N 
CY manifolds M vac (JV;X): 

X N 

f N (t;X) = f 1 (t;—). (7.8) 

on 

Here, we are again using the plethystic exponential relations. In fact, in cases 
in which the symmetric product is a complete intersection, f N will be finite and 
we can compare the computation of /jv using the formulas at the beginning of 
this section to independent derivations using the property that the manifold is 
a complete intersection. In cases in which the symmetric product is not a com- 
plete intersection we can still use the reasoning of §5.2. To count the number of 
generators and the number of defining relations for the symmetric product. 

Actually, eq. (7.7) is the reason why eq. (7.2) works. The far-LHS of the expression 
is the generating function for multi-trace at finite JV, in line with interpretation II stated 
above, while the far-RHS is the single-trace generating function at JV — >■ oo, in accord 
with interpretation I. Indeed, multi-trace operators with JV single-trace components in 
X is in one-to-one correspondence with single-trace operators in Sym N (X). Therefore, 
eq. (7.7) serves to bridge the two, whereby showing that the expansion coefficients Qn 
indeed count multi-trace operators at finite JV. Having delved into much abstraction, 
let us be concrete and now show how these proposals agree with known results. 

7.1 Example: The Complex Line C 

The simplest example, as was encountered in §6.1, is given by 

01(f) = /»(*; c ) = 7T37T = E r - ( 7 - 9 ) 
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This is the well known partition function of the half-BPS states in M = 4 SYM (given a 
particular choice of the supercharges). This partition function also counts the "extremal" 
BPS mesons in toric quivers (i.e. the mesons lying along an edge of the toric cone). In 
this case, it should be simple to check the multi-trace generating function is, as dictated 
by eq. (7.1) and eq. (7.2), given precisely by 

oo oo 1 

9(^t) = j: 9N (t)u-=]J-—-. (7.10) 

JV=0 m=0 V ' 

We note that is also the partition function of N bosonic one-dimensional harmonic 
oscillators. In other words, we are taking a quantum particle whose single particle states 
are precisely the integer points in the half-line Z> , and considering the placement of TV 
of such bosons. We can obtain g N for any value of N by Taylor expansion: 

N 1 

^(*) = n t— ^- ( ? - n ) 

71=1 

In fact, there is another way to see eq. (7.11). Indeed, we have the single-trace 
generating function explicitly: 

1 - t N+1 

f N (t; C) = l+t + t 2 + ...+t N = , (7.12) 

which encode the operators Tr(X*) for % — 0, . . . , N. We can take the PE of eq. (7.12) 
and using the multiplicative property eq. (5.5), arrive at eq. (7.11) directly. 
A few specific cases are at hand. For N — 1, 

gi{t) = ^— t = 1 + 1 + t 2 + . . . + t n + . . . , (7.13) 

corresponding to the operators 

I; TrpO; Tr(X) 2 ; Tr(X)"; .... 

For N = 2 we get 

9*(t) = -z ^ = l+t+2t 2 +2t 3 +3t 4 +3t 5 +4t 6 +. . .+(n+l)t 2n +(n+l)t 2n+1 +. . . , 

(1 - t)(l -t 2 ) 

(7.14) 

corresponding to the operators (we have dropped the Tr in the notation without ambi- 
guity): 

I; (X); (X) 2 ,(X 2 ); (X) 3 , (X)(X 2 ); (X) 4 , (X) 2 (X 2 ), (X 2 ) 2 ; . . . 
(X) 2n ,(X 2 )(X) 2 "- 2 ,(X 2 ) 2 (X) 2 "- 4 ,...,(X 2 )™- 1 ;... . 
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Indeed gz(t) is the Poincare series for C 2 /Z 2 where Z 2 acts by the exchange of the 
two coordinates (-21,-22) of C 2 . It has two generators, one of degree 1 corresponding to 
a := z\ + Z2 and another of degree 2, corresponding to b :— z\Zi- It is easy to see that 
the other invariant of degree 2 is represented in terms of these two, z\ + z\ = a 2 — 2b. 
Similarly, all other invariants of higher degree can be written in terms of these two. 
For X = 3, we have 

g 3 (t) = - r- 1 9W1 - = l+t+2t 2 +3t 3 +4t 4 +5t 5 +7t 6 +8t 7 +10t 8 +12t 9 +Ut w +O(t u ) , 

(7.15) 

corresponding to the operators 

I; (X); (X 2 ),(X) 2 ; (X 3 ), (X)(X 2 ), (X) 3 ; (X)(X 3 ), (X 2 ) 2 , (X) 2 (X 2 ), (X) 4 ; 
(X 2 )(X3), (X) 2 (X3), (X)(X 2 ) 2 , (X)3(X 2 ), (X)5; .... 

We see that indeed our generating function eq. (7.2) agrees with the explicit counting. 

To demonstrate the interplay between plethystics and symmetric products, we now 
calculate fx using eq. (7.8). We need to find /at(£; C) = fi(t; Now, if we expand 
an X-th order polynomial equation in one complex variable x, 

N N 

P N ( X ) =x N + J2 a i xN " = H( x - z i), ( 7 - 16 ) 
i=l j=l 

we find that the parameters Oj, i — 1 . . . N are symmetric functions of degree i for the 
coordinates Zj on C^. That is, are coordinates on .M vac (X; C) = C n /Sn- Further- 
more, there is precisely one generator of degree % for the ring of symmetric functions of 
the Zj for any % between and X. We can pick the generators to be the coordinates 
Oj. Any other symmetric function of degree i > N can be written in terms of the a^. 
Collecting this together we find fi for M vac (N; C) as in eq. (7.12), 

(T N 1 — i N+1 

/i(*;^) = n^- = /jv( * ;C) ' (7 - 17) 

consistent with our proposal from eq. (7.8) which implies, using the plethystic exponen- 
tial, eq. (7.7), thus supporting the proposal for multi-trace. 

7.2 Example: The Complex Plane C 2 

Next, we address a slightly more involved example, viz., C 2 . This case is quite simple 
as well and describes 1/4-BPS operators in U(N) M = 4 SYM. This also describes a 
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subsector of BPS operators in many toric quivers, namely the operators corresponding 
to points lying along a face of the toric cone (this face gives a toric subcone of the toric 
cone, a cone over a SUSY 3 cycle), when the SUSY 3-cycle has the topology of S 3 . Now, 
we have that 



hit; C 2 ) = 2t; 

oo 

9l (t; C 2 ) = PE\h{t- C 2 )] = -—- = j> + . (7.18) 
Formulae (5.6) and (7.1) takes the form 

oo ^ oo ^ k 

9( v > *) = II a _ utm)m+ i = ex pE kd-t^ (7 ' 19) 

m=0 V > k=0 V ' 

which looks deceptively similar to the generalized MacMahon function which is used as 
the partition function for the topological string on the conifold [66]. Using eq. (7.2), we 
get 

1 + 1 2 

g 2 (t) = — — — = 1 + 2t + 6t 2 + 10t 3 + 19t 4 + 28t 5 + 44t 6 + . . . . (7.20) 

(1 £)(! + £) 



We report the i?-charge 3 GIO's, corresponding to the term 10t 3 . We see that indeed 
there are 10 = 2 + 3 + 3 + 2 of them: 

(X 2 )(X), (X) 3 ; (X 2 )(F), (X)(XY), (X) 2 (F); [X ~Y] . 

Next, for i?-charge 4 GIO's, we see that there are indeed 19 = 3 + 4 + 5 + 4 + 3 them: 

(X 2 ) 2 ,(X 2 )(X) 2 ,(X) 4 ; 
(X 2 )(X)(Y), (X) 2 (XF), (X 2 )(XF), (X) 3 (F); 
(X 2 )(Y 2 ), (X 2 )(Y) 2 , (X) 2 (Y 2 ), (X)(XY)(Y), (X) 2 (Y) 2 ; 



The moduli space of vacua for this case is (C 2 ) 2 /Z 2 where the Z 2 acts as exchange of 
the two coordinates. It is a complete intersection and has 

h{t- (C 2 ) 2 /Z 2 ) = hit- C 2 ) = 1 + 2t + 3t 2 - t\ (7.21) 
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7.3 Example: The Conifold 

The single trace GIO's for the conifold are given, recalling eq. (4.5), by 

I __|_ 4- 00 

Jfc(*) = ?r f i yi = £( n + 1 ) 2 * B - ( 7 - 22 ) 

^ > n=0 

From formula (7.1) we get 

= l+t + 7t 2 + 3t 3 + 4t 4 = l+4t+19t 2 +52t 3 +134t 4 +280t 5 +554t 6 +. . . , (7.23) 
y ' (1 - t) 3 (1 — t 2 ) 3 

corresponding to the operators (again, we drop the Tr for brevity): 

I; (M M );.... (7.24) 

At i?-charge 2 we have 9 single-trace GIO's (cf. Fig. 4) and 10 double-trace GIO's, given 
explicitly by: 

(M 0) iM 0) i) (M ,iM 1)0 ) (M 1)0 M 1)0 ) 

(M 0) iM_ li0 ) (M 0) iM 0) _i)= (M 1)0 M_ li0 ) (M 1)0 M 0) _i) (7.25) 
(M_i )0 M_i )0 ) (M_i )0 M _i) (Mo_iMo_i) . 

and 

(M 0) i)(Mo,i) (M 0l i)(M 1)0 ) (M 1)0 )(M li0 ) 

(M 0) i)(M_ li0 ) (M ,i)(M 0i _i), (M 1)0 )(M_ 1)0 ) (M 1 , )(M ,_ 1 ) (7.26) 
(M_ li0 )(M_ li0 ) (M_ 1)0 )(M 0) _i) (M 0i _i)(M 0i _i) . 

We emphasize that all these GIO's have vanishing mesonic charge, we are not counting 
the BPS operators such as det(A). 

7.4 Refinement: Multicharges at Finite N 

As with the refinement of the charges discussed in§3.2.1, it is simple to generalize the 
arguments of the previous subsection to partition functions g± depending on more than 
one variables, arising for instance from CY cones with isometry U(l) 2 or U(l) 3 . 

Consider a toric CY cone whose integer points are described by the set C . The single 
particle states are described by the generating function 

9i(ti,h,t 3 ) = W*s- (7-27) 

n,m,r£C 
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Every point in C contributes once to gi, i.e., we are considering a quantum particle 
whose states are precisely the integer points in the toric cone. The multi-trace generating 
function <?jv(£i, £2, £3), in analogy with eq. (7.1), is given by 

1 / 00 k \ 

9M = J2 9 n(u)v n = n (i-^ S ) =exp ^T^'^'^ • (7 ' 28) 

The coefficients gN(U) can be interpreted as the multi-particle partition function of A r 
boson whose single particle states are given by the integer points of C. 

7.4.1 The Conifold Reloaded 

Recalling eq. (4.4), the generating function g 1 is given by: 

xy(l — q 2 ) , 1 l x 9 . 1 9 1 9 1 x y. 

' - l+g(a;+y+-+-)+g 2 (l+ — +:r 2 +— +y 2 +:r?/+ — +-+-)+. 



(1 - qy){q - x)(q - y) ' x y x 2 y 2 xy y x 

We can identify the charges (q, x, y) for following four operators 

M ,! = (1,1,0), M ,_i = (1,-1,0), M ,i = (1,0,1), M 0) _i = (1,0,-1) . 

Therefore, the generating function is, according to eq. (7.28), 

( ] J 1 1 

9{U ' q,X,V> ~ (1 " ") (1 " «/?*)(! " " " "J) (1 - uq 2 )(l - ^)(1 - 

1 

(1 _ Uq 2y2 ){l _ ^ )( 1 _ ^^(X _ _ ^ )( 1 _ ^S) • 

Now we try to apply above result. For N — 1 case we find the 

£1 = 1 + g ( - + :r + - + y ) + g 2 (l + ^ + a; 2 + ^ + — + - + - + ^ + y 2 ) + ... , 

y / V x y 2 xy y x J 

which has obvious correspondence with the variables Mjj. For = 2, we get the 
following expansion up to i?-charge two 

(\ 1 \ 9 / 2 9 2 2 2x 2y A 

#2 = 1+g - + + - + y )+q 3 + — + 2a; 2 + — + — + — + — + 2xy + 2y 2 . 
\x y J \ x 2 y 2 xy y x J 

Again, it is easy to find the mapping between the terms here and operators in eq. (7.25) 
and eq. (7.26). 
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7.5 Theories with only U(l) 2 symmetry 

We can also start from a theory with only U{1) 2 symmetry, whose Poincare series is 
given by 

giihM) = J2 a m,ntit2 • (7.29) 

n,m>0 

Using this we can find the generating function given by 

, / OO £ \ 

g^U) =^9N{t> N = II (1 . ytV w., = B P E^^) • ( 7 - 3 °) 

N n,m>0 V 1 2 ' \k=0 J 



8 Conclusions and Prospects 

In this paper we have considered the 1/2-BPS operators of generic superconformal quiver 
gauge theories, living on N D3-branes probing the tip of a Calabi-Yau (CY) cone. It 
is shown how to construct the explicit generating functions that count the scalar BPS 
operators. We have discussed in great detail various classes of CYs (orbifolds, toric 
varietes, del Pezzo's and complete intersections, even geometries for which the gauge 
theory is not yet known), providing a simple bridge (the "Plethystic Exponential") 
between the algebraic geometry of the CY and the generating functions of the BPS 
states. 

The plethystics directly relate three different generating functions: (1) the defining 
equations of the CY (syzygies) as well as the moduli space of vacua, (2) the single-trace 
operators and (3) the multi-trace operators. Beautiful structures thus emerge, exhibiting 
a rich inter-play between quiver gauge theories, algebraic geometry, combinatorics and 
analytic number theory. This intricate framework allows us to solve the 3 problems 
posed in the introduction, whereby realising our wish-list. 

There are a number of directions that could be pursued for future work. Let us dis- 
cuss some of them. We only considered the subset of operators with vanishing baryonic 
charges. For instance, for the conifold we did not include, in the counting, the operator 
det(A). It would be nice to find the partition functions including the baryonic charge, 
that may be compared to analogous computations on the string, AdS 5 x X 5 , side. 
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A possible continuation of our work could be to extend the study of chiral 1/2- 
BPS operators in M = 1 quivers to consider also 1/2-BPS operators with space-time 
angular momenta and 1/2-BPS fermionic operators. This would give partition functions 
depending on additional charges and would, for instance, enable a computation of the 
BPS index in quiver gauge theories, see [67, 68, 69, 70]. 

Another exension would be to consider 1/4-BPS operators, annihilated only by one 
supercharge. We remark that we are studying 1/2-BPS operators in M = 1 gauge 
theories annihilated by 2 out of the 4 supercharges Q. These are the analogues of 1/8- 
BPS operators (annihilated by 2 out of the 16 supercharges) in J\f = 4 SYM. It would 
be very interesting to extend the study to 1/4-BPS operators of quivers (annihilated 
by only 1 supercharge Q, analogous to 1/16-BPS ops in M = 4 SYM) 12 . One possible 
outcome, for instance, could be a comparison with entropy counting of the recently 
constructed AdS 5 SUSY black holes 13 . 

In M = 4 SYM an interesting problem is whether there is a change in the number of 
BPS operators changing the coupling. At zero coupling one expects more BPS states. 
Here, there is a precisely analogous question. It was shown in [28] that, in the moduli 
space of SCFTs corresponding to a given quiver gauge theory, there is a special point 
with enhanced chiral ring 14 . It was observed in [28] that at this special point, the 
growth of the number of single-trace, N = oo, BPS mesons is exponential instead of 
that quadratic (as is the case on generic points of the moduli space of SCFTs). It would 
be interesting to study further this mechanism, that could lead at finite temperature to 
phase transitions. 

On the gravity/string side of AdS/CFT, we should also find the same partition 
functions. For single-trace operators the result is well-known. The interesting case is 
multi-trace at finite N. One way to reproduce the g^ should be counting Giant Gravitons 
(GG's) in AdS 5 x X 5 . There are studies of GG's in AdS 5 x S 5 and AdS 5 x T 1,1 [71]. In 
the case of S 5 , one considers the classical moduli space of GG's and is led to study N 

12 One single-trace example of such operators is given by Tr(0 K) 1 where O is a scalar BPS operator 
and K is the scalar SUSY partner of a conserved current: Q a O — 0, and Q 2 K = Q 2 K = 0, so Tr(0 K) 
is annihilated only by Q 2 . 

13 Notice that these BPS black holes are constructed in 5-dimensional gauged supergravity with U (l) fc 

gauge group, so can in principle be uplifted to various AdS$ x X 5 solutions. 

14 For N = 4 and orbifolds thereof this special point is the free theory, for the conifold it corresponds 

to having vanishing superpotential. For generic quivers we have only one term of the superpotential 

vanishing. 
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classical non- interacting particles, whose single particle phase space is C 3 . Quantizing 
the multi particle phase space one gets iV bosons whose single particle states are the 
integer points of the toric cone of C 3 . The partition function is precisely the finite N 
partition function of 1/8-BPS ops in Af = 4 SYM. In generic toric quivers, for instance, 
one should find the result of the end of §7: 



A different approach is [72], where they consider "dual GG's", i.e., D3-brane wrapping 
an S* 3 inside AdS^, and moving along a trajectory in X 5 . For generic Sasaki-Einstein 
manifolds X 5 , we conjecture that these states are BPS if and only if the trajectory is a 
BPS geodesic. We already know that single-trace BPS mesons are the quantization of 
such BPS geodesies (see [28]), so the final result should be a simple outcome of the use 
of the Plethystic Exponential on the gravity side. 
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